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SYNOPSIS 


This thesis entitled "OPTICAL PROPERTIES OF RANDOM 
SEMICONDUCTING ALLOYS" is submitted by Biplab Ganguli to the 
Department of Physics, Indian Institute of Technology, Kanpur in 
partial fulfilment of the requirements of the Ph. D. degree. 

The problem of optical responses likeoptical 
conductivity, dielectric function etc. of disordered alloys is a 
matter of interest for a long time now both for theoretical and 
experimental physicists. The optical response of ordered systems is 
well understood on the basis of the formulation based on the Bloch’s 
Theorem and band structure ideas. The electron photon interaction 
matrix element factor in the expression for optical conductivity is 
assumed to be slowly varying function of energy. It is found that 
for a large class of materials the electron -phot on interaction 
matrix element is a slowly varyig function of energy over a large 
energy range and fairly independent of randomness. However, the 
need to go beyond these assumptions for disordered alloys was felt 
for a long time. 

In Chapter II, we discuss optical response in general and 
find a basic formula for optical conductivity. In particular we set 
up expressions for the optical conductivity for ordered systems and 
disordered alloys for two cases, one with the assumption that 
electron -photon interaction matrix element is independent of energy 
and randomness and the other without this assumption. 



( viii) 

In chap'ter III, we discuss a self —cons i stent, coherent 
potential approximation CCPAD and cluster coherent potential 
approximation CCCPAD formalisms to obtain a configurationally 
averaged Green function. The cluster CPA formalism is based on the 
AUGMENTED SPACE METHOD introduced by Mookerjee C1973D. The 

Augmented space method is by now a well established and powerful 
formalism for calculating various average properties of disordered 
systems. In CPA the correlated scattering from two or more sites 
are ignored. In the self-consistent cluster generalization CCCPAD 
of CPA statistical clustering effects and off-diagonal disorder over 
finite cluster sizes are successfully taken into account. Beside, 
it preserves the CHerglotzD analytic properties of the Green 
function which are essential for obtaining physically meaningful 
approximate results. 

In this chapter the Augmented space method is used to 
obtain an expression for configurationaly averaged optical 
conductivity of random semiconducting binary alloys without making 
any assumption about the electron-photon interaction matrix element. 
Thus a genral formulation for optical conductivity of random 
semiconducting binary alloys is obtained. This genralization is the 
principal contribution of the thesis. The formulation is applied to 
a model case. We use CPA and CCPA to calculate density of states 
and the various optical response functions. The effects of various 
assumptions are also discussed within the model. 

In chapter IV we make applications of the above formalism 
to realistic systems. The systems chosen are GaAs. Sb and 
Ga^_^In^Sb, In order to compare with the older formalisms we 



ix 

calculate various optical response functions for the two cases. 
Firstly using the above formalism and secondly making assumptior 
that the electron -phot on interaction matrix element is randomness 
independent and taken to be the averaged value. We also calculate 
the square of the electron-photon interaction matrix element, and the 
product of this and the valence band density of states as a function 
of energy to accurately interpret the energy resolved photo-emission 
experiment. The calculations are done for two concentrations, i . e. , 
X = . 5 and x = . 2. 

In this chapter we also discuss molecular orbital basis 
Cbond basisD which is used to generate Recussion coefficients for 
the calculation of valence band density of states. TTie Recussion 
method Introduced by Haydock C1972D enables us to work on such 
realistic three-dimensional lattices. We discuss Recursion method 
in this chapter. The Recursion method is a well-established means 
of solving the Schrddinger equation for a system in which local 
interactions dominate. The local density of states and related 
quantities are calculated from the continued fraction expansion 
generated through the application of the Recursion method. We also 
discuss vaTious terminator schemes to terminate the infinite 
continued fraction expansion. The Recursion method together with 
CPA is used to calculate the Green function, the density of states 
and the various optical response functions. 

In chapter V, the achievements and shortcomings of the 
present work are discussed. Some of the problems in this field are 
also suggested. 
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CHAPTER I 


INTRODUCTION 

The problem of optical response in disordered alloys has 
received the attention of solid state physicists with the extensive 
experimental works in that field for the last two decades CNilsson, 
1974 ; Zavetova & Velicky 1976D » The theoretical understanding of 

optical response of periodic systems is well understood with the 
application of Bloch’s theorem and band structure methods. 

Optical conductivity is related to the product of the 
Joint density of states and the electron -phot on interaction matrix 
element Coptical transition matrixD. For a large class of ordered 
systems, it is found that the optical transition matrix element is a 
very slowly varying function of energy and, thus, can fc>e thought of 
approximately as a constant. Therefore the optical conductivity is 
directly proportional to the product of the frequency of the 
incident photon and the integrated Joint density of states. 

The theoretical understanding of optical response of 
disordered alloys is not so satisfactory. The reason for this 
inadequacy is an additional problem related to the study of 
disordered systems : configurational averaging. In the case of: 
ordered systems, once the Hamiltonian model is set up, the solution 
of quantum mechanical problems is essentially the solution of 
Sfchrddlnger equation. However, with the introduction of disorder,! 
we face a new problem, i.e. the lack of sufficient Information about: 
the sample in the Hamiltonian. litis itself involves random parameters 



and until information about the distribution of these parameters is 
supplied, the solution remains incomplete. One of the central 
problem in the study of disordered systems is configuration 
averagi ng. 


Before going on to say anything about how we tackle the 
problem of configuration averaging, let us first understand what the 
term coTiftgva:-a.tion implies, why configuration averaging is needed 
and what is it that we should average. In any substitutionally 
disorder system, the lattice sites are occupied by random potentials 
: random in space as in quenched disordered solids or random in time 
as in thermally disordered systems, Cfor example, electrons in 
contact with a phonon bath? . There are examples where both spatial 
and temporal disorders are involved, for example dirty alloys at 
high temperatures CMooij systems!). A simple example of a spatially 
disorder system is a random binary alloy. Cu-Ni is a good 
candidate. At T = the randomness is quenched and spatial. The 

potentials are characterized by random parameters. For example in a 
random binary alloy the potential, VCr^D = * VgCl-N^D, is 

characterized by the random, occupation number C = O if i‘*': 

site is occupied by an A atom and = 1 if i*"^ site is occupied by! 
a B atonD. If short ranged order is absent, then the probability! 
that =» O or 1 is directly proportional to the concentrations of| 
the A or B species. A particular realization of these parameters,; 

j 

either in a given sample or at an Instant of time is what we refer 

i 

! 

to as a configuration of the system. | 

i 

! 

! 

Let us now understand why it is necessary to take| 

i 

configuration averaging. We shall start with the above example of &' 


random binary alloy. A configuration is determined by a particular 
assignment of O’s and I’s to the sequence <N^>. Any one sample 
therefore corresponds to a particular configuration. For a sample 
containing N atoms there are 2*^ configurations. This is a very large 
number for a macroscopic system. An experimentalist at the outset 
has no knowledge of the particular configuration of his sample. Nor 
is such a knowledge of much use. Evidently, every different sample 
of the alloy will give experimental information having different 
microstructure according to the particular configuration it has. 
The experimentalist, on the other hand, is primarily interested in 
the overall statistical trends of his results which he interprets as 
the physical propertiGS of the random alloy. A theorist too can 
generate a number of results depending upon the particular 
configuration he chooses for his random parameters. This again is 
not very useful. He could try to derive properties of the most 
probable configuration and propose that these properties are highly 
likely. But, again, the number of different configurations as 
probable as the most likely and those with comparable likelihoods of 
occurrence, are intractably numerous. Therefore, we have to evolve 
a statistical description of the system. The most basic of these 
descriptions is the averaging over all possible configurations. 

It is clear from the above why configuration averaging is 
required. It is also clear that it is the physically measured 
quantities which must be configurationally averaged. It was pointed 
out by Anderson in 1972 that the type of statistical description 
essentially depends on what we are trying to describe. He warned 
that only experimentally observable physical quantities related to 



the response of a. physical system to a probe must be averaged. Thi: 
warning was necessary to explain the source of controversies thai 
arose over non-physical results obtained because directly 
non -measurable quantities like self -energies, wave -functions etc. 
were being averaged. Whereas, it is quantities like the density of 
states, and response functions like optical conductivity or 
diffusion probabilities which should be averaged. Since the 
imaginary part of the Green function is directly related to the 
density of states. Green function techniques are favoured for the 
study of disorder systems rather than the Schrddinger equation 
appr oach. 


Configuration averaging is a valid procedure only when the 
probability distribution of the property under consideration is 
sufficiently well behaved so that the average dominates over the 
higher moments. Interesting examples are the intensity of starlight 
transmitted through layered media with randomnly varying refractive 
indices studied by Chandrasekhar C1960D and the problem of the 
resistance of disordered chains. Kumar and co-workers Cl 985, 1986;) 
have shown that the variance of the resistance of disordered chain 
diverges much faster than the average as the length of the chain 
increases. It is, therefore, not meaningful to talk in terms of the 
average resitance of a very long disordered chain. 

Before we begin the discussion of how configuration 
averaging may be carried out, it is worthwhile to say something 
about spatial ergodiclty. For a large enough systems, in the limit 
of infinite size, all possible environments are achieved even in a 
single sample. The figure Cl. 15 illustrates this idea. 


A global probe samples all possible envlronmen'ts in the 
sample weighted with appropriate probabilities. A global property 
is thus automatically configuration averaged, even if it refers to a 
single sample. This sampling of the entire configuration by a 
spatially global probe is analogous to the sampling in statistical 
mechanics of all configuration space by a dynamically evolving 
system. Hence the name -.spatial srgoclicity. Such global properties 
need not be configuration averaged all over again. This idea is 
implicitly behind a large body of the 'Targe cluster" calculations 
involving either isolated or immersed clusters. One example is the 
application of the Recursion method to large disordered clusters 
CHaydock, 19871). 

In general , the exact calculation of configuration 
averages is an intractable problem. In §3.2 in chapter III, we shall 
discuss in detail how to configuration average any physical quantity 
in a simple tight -binding model for a substltutionally disordered 
alloy. We shall discuss the configuration averaging of optical 
conductivity and related response functions. 

Various mean field approaches, like the self-consistent 
coherent potential approximation CCPADCSoven, 19673, provide 
powerful and physically meaningful methods for the study of 
configurationally averaged one and two particle Green functions for 
disordered systems which are directly related to the averaged 
density of states and response functions. The earlier works on 
optical properties of disordered systems CZ^v^tova and Velicky 
C19763 ; Cody 1984; Pickett et.al. 1983; Zdetsln et. al . 1985, etc. 3 



used the mean field and other approaches without taking into account 
the vertex correction arising out of disorder dependence of the 
optical matrix element . That is, while doing the averaging, the 
average of the product of density of states and optical transition 
matrix element was replaced by the product of their averages. 
Pickett and co-workers C1983D used the virtual crystal approximation 
CVCAD to average the optical transition matrix element separately. 
They, then, took this factor as constant and finally used coherent 
potential approximation to get averaged density of states. 

This procedure may suffice for large classes of systems, 
but for more accurate studies the full problem has to be tackled. 
Particularly, in the strong scattering cases, the vertex correction 
may be of importance and the optical transition matrix element may 
not be insensitive to energy variations. 

The need was felt, therefore, to include the vertex 
corrections. Mookerjee C1973a, b, 1975 a,b,c5 introduced a 
formalism called the Augmented space formalism CASFD to calculate 
configuration average of a general function of many random 
variables. The formalism has been quite successful in making the 
generalization of single site CPAto Cluster Coherent Potential 
Approximation CCCPAD and also to include the off-diagonal disorder 
in a self-consistent manner C Kumar et.al. 19823, Nickel and 
Krummhansl Cl 971 3 tried to make such a generalization using a method 
based on the corrected cumulant expansion to find electronic density 
of states of a one dimensional model of an alloy. Butler Cl 972, 
19733 calculated the real and imaginary parts of the Green function 
using a self-consistent cluster method. These generalizations posed 



analytic problems, such as multivalued, discontinuous and negative 
density of states in certain energy regions. 

Haydock e t . a I . Cl 972D , Mook er J ee Cl 973a , bD , Mtil 1 er -Har tmann 
C1973D and Mills and Ratanavararaksa Cl 9785 pointed out that for any 
real potential function, the Green function must obey certain 
mathematical properties called the Her-glotz properties. A complex 
function fCz5 is defined as Herglotz if 
Ci5 ImfCzD < O for Imz < 0 

Cii5 Singularities of fCz5 lie on the real axis, and 

Ciii5 fCz5— > ^ as E — > oo Cwhere z = E + iO'^5 

In order to get physical results, Herglotzici ty has to be 
retained in any approximation to the Green function. It has to be 
borne in mind that in the development of a theory which Includes 
scattering from clusters, short ranged order, off-diagonal disorder 
or the effects of positional disorder, Herglotz property must be 
inherently retained. 

The Augmented Space Formalism allows us to systamatically 
introduce approximations which retain such properties of Green 
function. The main aim of the present work is to extend the 
Augmented ^ace Formalism to obtain configurationally averaged 
optical properties including the vertex correction. 

Chapter II deals with the detailed discussion on general 
optical properties of solids, response functions and their relations 
to one another. In chapter III we Introduce the Augmented ^ace ' 
Formalism with its mathematical formulation. It is shown how we can ; 
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generate CPA and CCPA within the formalism. In §3. 6 a general 
expression for the optical conductivity of random semiconducting 
alloys is formulated using the Augmented space formalism In §3.7 we 
apply the formalism to a model case. Chapter IV deals with the 
application of the formalism to the ternary systems GaAs^ Sb and 
Ga^_^In^Sb for different concentrations c. In this chapter we also 
discuss molecular orbital CMO:> or bond basis and Recursion method 
and its different terminator schemes which enable us to 
accurately calculate the constituent Green functions. These are 
used as essential inputs in our theory. Side by side, for 

comparison, we also calculate the various response functions without 
the vertex correction. 



Representation of different environments 

within a sample. For n > 00 , all 

environments are realized. 


Figure 1 . 1 




CHAPTER II 


THEORY OF OPTICAL PROPERTIES 

2. 1 INTRODUCTION : 

The interaction between electromagnetic radiation and 
solids causes absorption of the radiation through microscopic 
excitations. The excitations can usually be described with 

quasi par tides like electrons, phonons, plasmons, magnons etc. For 
metals, the dominant absorption is caused by the electrons. In the 
infrared CIRD region CO. 1 eVD the conduction or "free" electrons 
contribute to a strong absorption, the magnitude of which decreases 
smoothly with increasing photon energy. In the visible region C2 to 
3 e\0 it is usually comparatively small C Nilsson 19743. The 

microscopic process can be viewed as a electron-photon scattering 
process. The scattering occurs to a conduction or unoccupied band. 
Such excitations are called Intraband transitions. This process is 
shown in the figure C2. 1.13. The theory for such an absorption 
process was first given by Drude and the model is known as the Drude 
model CAnimalu, 19783. In his model, Drude used a classical 
approach. The basic assumption was that the optical conductivity 
and the dielectric function could be determined by considering the 
motion of quasi -free electrons under the influence of the 
oscillating electric field vector of the electromagnetic wave. 


The optical spectrum of a real crystal of a metal , a 
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semiconductor, an insulator or a disordered system usually has mor 
structures in the visible region to far UV region CIO eV5 than thos 
arising from E)rude absorption and the occurence of plasm, 
oscillations of free electrons accelerated by the electric field oi 
the light wave. These structures are seen in the form of peaks, 
edges, etc. , in the absorption versus frequency curve. SucJ 
structures arise due to the interaction of electrons only with the 
radiation. The electron is then directly excited from one 

electronic band to another. The process is called Interbanc 
transition. This priDccess is shown in the figure C2. 1.2Z>. Ir 
materials like Insulators, semiconductors, or alloys of these, this 
is the only kind of absorption that takes place, because in sue! 
materials the valence band is completely filled and excitation car 
take place only to the completely unfilled conduction band. Drude 
absorption, therefore, is not possible. To understand these 
processes, we need to go beyond the classical approach and develop a 
quantum transport theory. In such a theory the electron -phot or 
interaction is considered to be the exciting part of the 
Hamiltonian and a time dependent perturbation formalism is appllec 
to obtain expression for the optical response functions of the 
system. Such a linear response theory was first suggested by Kube 

C19573. ! 

! 

2.2 OPTICAL RESPONSE FUNCTIONS ; i 

In this section we shall study the linear response of 

1 
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electrons in a solid to external electromagnetic radiation, from the 
IR to the UV region CO. 1 - 10 e\0 . In particular we shall describe 
the dielectric function sCuO , optical conductivity crCox) and 
permeability pCoO as the main properties usually studied 
experimentally. Here to is the angular frequency of the external 
electromagnetic radiation. In the optical region CIR to U\0 we 
shall choose non-magnetlc systems where 2 = 1 . The dielectric 

function sCoD , conductivity function oCca) and permeability function 
fjCoD are defined through the response disturbance relations 


DC 0)3 

== sCoD 

EC 0)3 

C2. 2. 13 

JCo>3 

- cK 0)3 EC 0)3 

C2. 2. 23 

HCo)3 

_ 1 

BCo)3 

C2. 2. 33 

^iCo)3 


respectively, here D, E, J, H and B are the displacement field, the 
electric field, the current density, the magnetic field and the 
magnetic induction vector respectively. A related response is the 
refractive index function nCuO which is related to sCo£> through : 

nCciO = C2. 2. 43 

Another related response is the absorption function aCttO.: 
It is defined as the fraction of energy absorbed! n passing through 
the unit thickness of the material. 

Electromagnetic waves are the solutions of Maxwell's 
equations. For linear and homogeneous media these are written as 
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follows CIn CGS units!) 


Curl H = L ^ + lEj 
c dt c 


ca. a. 55 


7. D = O 


ca. a. 65 


Curl E = 


1 ^ 
c ^ 


ca. a. 75 


7. H = O 


ca. a. 85 


Eliminating the magnetic field vector between these Maxwell ’s 
equations for the electric field vector : 


7*E = 


c d E . Ana HE 


c* i>t* 


2 dt 


ca. a. 95 


This represents a wave propagating in a solid with dissipation Csee 
figure Ca.a. 155. If we choose the electric field of the type 


E = E^ exp <iC«.r - wt5> 


ca. a. 105 


then our wave equation requires 


4no'i oj 


i . e. , 


le = 


r . . 4 tio't 

[« * i -TT ] 


ca. a. 115 


In general the propagation constant a comes out to be a 
complex number. In free space we would simply have 


ca. a. 135 
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Figure 2.2.1 Damping of Electromagnetic Waves in Solid 
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In our medium the wave is travelling with the modified velocity v. 



where n is the refractive index of the medium. Therefore 


X = ~ n Ca. 2. 13D 

c 

Comparing the relations C2.2. IID and C2. 2. 13D we find that the 

^Sijp 

refractive index nCwD comes out to be a complex function and is 
gi ven by 


nCtoD = + i C2. 2. 143 

Comparing the relations CE. 2.4D and C2. 2.14]>, a complex dielectric 
function is defined as : 

sCoD = + i - C2.2. 15D 

to 

The real part of the complex dielectric function eCoD is the 
measured dielectric response of the medium and the imaginary part is 
related to the optical conductivity. 

We may write nC<o> as 

nCto3 = nCtoZ) + in C<o5 C2. 2. 16D 
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where nCwD Is the measured refractive index. Then 


aee 




nCoD + in CoD 
2 


.A. 

e 


where is *the direction of polarization of electromagnetic wave 
Therefore the wave equation C2. 2. lOD becomes 


E^exp ^ ^ 


ne . r 

r 


^ ri coe . r^ 

- *■]] (- 


C2. 2. 171) 


We see that the velocity of propagation in the medium is reduced to 

tZnn 

— I 

per wavelength. The damping of the wave is, of course, associated 
with the absorption of electromagnetic energy. To calculate this we 
should use Maxwell’s equations to find the current associated with 
the equation C2. 2. 17D. This is the right hand side of the equation 
C2. 2. 5:> 


or , 




iW ~2 = 
— n E 
c 


C2. 2. 18D 


where is a complex current density. The rate of production of 
Joule heat is the real part of 


J^.E = - if |E|* 


ce. a. 1S3 
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Thus t-he absorpt-ion coefficient aCuO is given by 


or , 


aCcoD 


ReCJ .ED 
c 

n|E|" 


2n oj 
2 

C 


aCoiD 


4n<yCc>)D 

nCcoD 


C2. 2. 20D 


The optical absorption coefficient is perhaps most 
susceptible to simple experimental determination. Experimentally, 
it is measured as follows: when a beam of light of intensity I 

o 

goes normally through a slab of a medium of thickness x, it 
attenuates in accordance with the exponential law, I = I^expC-ooO. 
Consequently, a can be obtained simply by measuring of* ^he 

intensities impinging and emerging from the sample. 

It is seen that the whole game of optical response can be 
described by a single complex response function, viz -the complex 
dielectric function ^CcoD. The real and the imaginary/ parts are not 
independent. The analytic properties of cCwD lead to the 

Kramers-Krdni g relation: 




_ " Ca’-uOc CojO 

^ r i dW 

n J ,22 

O CO -<o 


C2. 2. 21D 


and optical 


sum rules 


00 


J to/p CcoDdoj 

2 



C2. 2. 22i 
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J , ji j[ 2 

o> Im s CwDdco = =Trco 

d, p 

o ^ 

AiT 

where s = s + and is the plasma frequency. 

We have already seen, the dielectric function relates the 
displacement current vector D and electric field vector E. But in 
rapidly varying fields, such as the fields associated with 
electromagnetic radiation, the field strengths involved are fairly 
small. Now any variable field can be reduced by means of a Fourier 
expansion to a set of monochromatic components with a time 
dependence e For such a field the relation between DCtD and 

ECtD can always be taken to be linear. The most general linear 
relation between DCtD at all previous instants can be written in the 
integral form 


_ _ ® _ 

DCO = ECO , f fCOECt-rtdr CE. B. 833 

o 

where fCzD is a function of time and of the properties of the 
medium. Taking Fourier transform of the equation C2. 2.S3Z) and 
comparing with the equation C2. 2. 15, we get 


£:C w5 =1 


00 . 

+ J'fCTDe^^’^dT 
o 


C2. 2. 245 


Thus for periodic fields the complex dielectric function! 

i 

I 

may be regarded as a function of the field frequency as well as of ! 
the properties of the medium. The variation of sCoO on the 
frequency is called the dispersion relation for the medium. The I 
function sCaD is in general complex for conducting media. From! 


the equation CS. 2. 24;> we see that 


= «*Co)D C2. 2. 25D 

Separating the real and imaginary parts, we have 

£C-(oZ) = £CcoZ> CS.S.SfSy 

and £ C -toZ> = -£ C w!> 

2 2 

Thus £CoD is an even function of the frequency and 
is an odd function. 

2.3 INTRABAND AND INTERBAND TRANSITICa^S : 

In section C2. 13 we have already discussed about the two 
kind of transition processes : intraband and interband, and seen the 

j 

distinction between these two. In this section we shall talk more 
about Interband transition. From figures C 2. 1.13 and C2. 1.23, it is 
clear that intraband transition is possible only in a metal whereas: 
interband transition is possible for all materials. i 

i 

Inter band transition takes place in the range of; 
frequency o>^ < <»> < where : 

to is the threshold value for inter band transi tion-lies in the: 
o I 

I 

visible region, and | 

j 

to is the plasma frequency-lies in the far UV region. j 

P j 

i 

The quantum transport theory for Interband transition can| 
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be developed for two different cases. First, at the temperature o'^K 
and the other, at some finite temperature. The nature of inter band 
electronic transition, in these two cases, is very different from 
each other. Since at o'^K there is negligible number of phonons 
present, the electrons absorb only photon energy. 



The 

momentum conservation requires 





hk^ 

= hk^ + hQ 


C2. 

3. 13 

where k^ 

and 

k^ are the final and initial wave 

vector 

of 

the 

el ectron 

and 

Q is the wave vector of the photon. 

But 

we 

are 


constraining ourselves only in the region of IR to UV light for 
which wavelength is much large compared to the lattice constant. 
Therefore, Q is much smaller than the k’s. Thus the momentum of the 
photon may be neglected and finally we get 

^f=^i 

or k_ = k . . C2. 3. 33 

f i 

This shows that transition takes place vertically upward in the 
reduced Brillouin zone scheme as shown in figure C2.3.13. This kind 
of interband transition is called direct transition. The energy is 
also conserved during such a transition 

E Ck 3 - E.Ck.3 = hw 

f f t I 


C2. 3. 43 
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The other kind of inter band transition is the Indirect 
:,ransltlon. This kind of transition process is observed in the 
ligher temperature region. At finite temperatures there are 

sufficient phonons present. Therefore electron absorbs phonon 
snergy as well as photon energy. Consequently, in this case, the 
::onservation of momentum gives 

hk^ = hk^+ hQ + hq C2.3. 5:5 

where q is the phonon wave vector. 

Since hQ is O 

therefore ^f ~ ^i ^ 

Hence k^ k^ C2.3. 6:> 

This process is shown in figure C2. 3. 25. This means the transition 
is not vertical in the reduced zone scheme. 


2.4 BASIC FORMALISM FOR OPTICAL CONDUCTIVITY : 


The fundamental expression for optical conductivity can be 
derived simply from first order time dependent perturbation theory. 
The crystal Hamiltonian is treated as the unperturbed Hamiltonian and 
the eigenfunctions of this Kami 1 toni 'lan as the unperturbed 
eigenfunctions. The one electron Schrbdinger equation for such 
unperturbed system is then given by 





83 



ca. 4. 13 


where H : crystal Hamiltonian 

VCFd : crystal potential 

V'*^CF>: eigen functions of H 

k 

m : the effective mass of the electron in the solid 


When electromagnetic radiation falls on the material, the photons 
interact with the electrons in the solid. This electromagnetic wave 
carries a usually weak electromagnetic field. The electromagnetic 
radiation adds an extra potential to the crystal Hamiltonian during 
interaction. This weak, time dependent potential is treated as the 
perturbation. The full Hamiltonian is then given by 

H = H + H C2. 4. aD 

— — O —4 

In general the Hamiltonian associated with an electron in a 
electromagnetic field is 

H' = 5 ^ C2. 4.3D 

where p = ^ V is the momentum operator 

ACr,t5 ; Vector potential associated with the 
electromagnetic field 
t3 : Scalar potential associated with the 
electromagnetic field. 

In the absence of free charges we may choose 
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<^r,t5 «= O 

Thus H* reduces to 


H’ 


2m»* 


CP 


e 

c 


ACr.t:^* 


C2. 4. 43 


C2. 4. 53 


or 


H’ 


Cp.p + — A. A - — p. A 
2m« ^ 2 c ^ 

c 


A. p3 


C2. 4. 63 


If ACr,t3 is a weak potential, to the first approximation the second 
term in the equation C2.4.63 may be neglected, leading to 


H": 


P- P _ 


2m^ 


® <p. A + A. p3 


C2. 4. 73 


Therefore H ia now given by 


H = - 


2m»tc 


Cp. A + A. p3 


C2. 4. 83 


Absence of free charges leads to 


V. A= O 


C2. 4. 93 


Since 


[A, p3 == ihV. A, therefore, p. A = A. p 


Putting back A. p in place of p. A in the equation C2.4.83, we get 


«i = 


m»tc 


A. p 


C2. 4. 103: 


This extra term is responsible for excitations. Electrons in an 
unperturbed state v'^Cr3 at a time t = O will subsequently have a 
wave function which can be linearly decomposed into the 
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complete set of eigenstates. <v/°CrI>> of H^. Here v/^cFl) denotes the 
state before the transition or initial state and v/^cFd is the state 
after the transition, i . e. , final state. 

The transition rate between an occupied state and 

unoccupied state may be calculated by solving the equation for 
vCr.tD by the standard first order time dependent perturbation 
theory. For a monochromatic vector potential. 

A = A exp ClCwt:>] C2. 4.1i:> 

o 

the rate of transition accompanied by the absorption of a photon of 
frequency co from the radiation field is 

Q = ^ l<V^-|H I ^>1^ <5CE -E.-hwD C2.4.12:) 

k k. 
f x. 

The energy absorbed per unit volume per unit time as 

result of this transition is Q , where O is the volume of 

k-.k. 

f t 

the sample. The total rate of energy loss is obtained by summing 
this over the initial and final states of the electrons in the 
system. Introducing a factor f,Cl-f,3 for this purpose In equation 
C2. 4.1E>, where f^ is the probability of occupation of the state 
V/^cFd and f^ is the probability of occupation of the state v<j.cF, t3 
we find that the total rate of energy loss is 

S ?? «r-r- 
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- 1 rs f.ci-f,3 t>o> §2 |<V',|H |v'.>|' «<:E^-E.-h<^ 

t f 


2n(o 


^ ZZ f .Cl-fp |<V',K— §-A. p5 |v^.> I* 6CE-E.-fi6>D 
fl . , V f • • m*«c I I f ^ 

I { 


CO ZZ f.Cl-fp l<V',|A. piv'. > 1"" <5CE -E.-hcoD. 

m^c^O it ^ ^ ^ 


But this is also equal to the joule heat loss due to the 
mobile electrons, i.e. J . E. But J. E = c/E. E, Cfor linear and 
i sot r opi c medi uirO . 

Since E = — A , therefore 

c 

T.E = <y CiOk>/c:>* A. A 


or 


J. E = - O' Cco/cD^ A. A 


C2. 4. 13D 


Equating these two, we find 


T — _ Sire 

a Cco/cj A. A = 

2 2 
c 


— j<Vf |A. p|v^.> I V.Cl-fp 6CE^-E.-h<oD 


If A = Ae^, where e^ is the polarization direction of the 
electric field vector E, the above equation becomes 


oCcoD 


3rre 1_ 
m»t*n t f 


EE l<Wf|e^-P|V'.> 1* <5CEj.-E-hcoD C2.4.1 


43 


This is the basic formula for the optical conductivity of a solid 
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due to inter band transition. This can be used for any kind of 
material like metals, semiconductors, or insulators and in a 
slightly modified form ideally suited for disordered alloy systems. 


2.5 ANALYSIS OF THE CONDUCTIVITy FORMULA: 


Looking at the basic formula for the optical conductivity 
given by the equation C2. 4.14D, we can have some idea about the 
physical behaviour of the optical conductivity for both ordered and 
disordered systems. 

CHEERED SYSTEMS 


There are two important factors in the expression C2.4. 14D 
for the optical conductivity, oCoiO : 

CID Conservation of energy factor : 6C E^ -E^^ -hwD 

C2]) Square of the optical matrix element : Kw, P | • 

Let us first denote the optical matrix element as 

X,tf 

We shall now show that the first factor is related to the Joint 
density of states. 

The one-electron density of states C number of allowed 
energy levels per unit energy range and F>er unit volumeD is given by 

nCE> = 2n“* E <5CE-E > CE. 5. ID 

n 

n 
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where O is the volume of the specimen. The factor 2 arises because 
of summation over the two spin orientations. From the properties of 
Dirac <5- function: 

J 6CE-Ej,5<5CE^-ED dE = 6CE^-E^D C2. 5. 2D 


We have 


6CE -E.-hcoD 

f t 


_max . 

r ^ 

= J <5CE-E p6CE-E.-ha>D dE 

Oft 


C2. 5. 3D 


where is “the upp^er limit, or conduction band. Therefore 

using equation C2.5.3D we have 


E E <5CE -E.-hooD 

T . ft 

f t 



or 


E E SCE-E^-hoD 

. ft 



max 

f 


< E 6CE-E D > < E 6CE-"E.-"ho>D > dE 
f i '■ 


C2. 5. 4D 


Using the definition C2. S. ID, we can write the right hand side of 
the equation C2. 5. 4D in terms of the conduction band density of 
states, n^CED and valence band density of states, n^CED : 


—max. 

Cl r 

r x: 6CE~E.-ho£> = J n CED n.CE-hwD dE. 

~~ fi 40f t 

f T. 


But nCEDn.CE-hwD = J,.CE.«D, where J.CE,wD is the Joint 
f I fv fl 


density of states. Therefore 
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—max 

E E < 5 CE -E.-ftoO = — J J, .CE,<.j 3 dE 
“T fi *0 ft 

f t 


C2. 5. 55 


For an ordered system, where crystal potential is periodic 
and the . wave-function is described by Bloch's functions, the 
optical matrix element factor is usually a very slowly varying 
function of energy CZavetova & Velicky, 19765. Therefore in the 
first approximation t>e replaced by its k-average, and 

takenout of the integral in the equation C2. 4. 145. Within this 
approximation this equation becomes 


CfCox) 


2fTe^ 1 _ 
m»t*0 ^ 



E E 


6 CE -E.-hw 5 

f I 


C2. 5. 65 


For semiconductors, at T = 0°K ,f. =1 and f, = O. Therefore in 

u f 

this case equation C2. 5. 65 becomes 


aCiiD 


Sne^ 



E E < 5 CE -E.-hoi 5 

r 7 f I 
I f 


C2. 5. 75 


Using the relation C2. 5. 55 , the equation C2.5. 75 reduces to 


-max 


oCa >5 


2ne 

m*f*n 


L 


rZ Q 


H f y 




J,.CE,co 5 dco 


C2. 5.85 


Therefore for a fixed value of to, we have 


to o<to 5 oc J., CE,to 5 dE 

Lf 


C2. 5. 95 
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This shows that the optical conductivity for any ordered system is 
directly proportional to the integrated Joint density of states. 

In an energy resolved photoemission experiment one 
essentially studies the energy distribution, NCE,CjOiD, of photoemitted 
electrons for a fixed value of frequency of radiation as a function 
of E. This is done by introducing a filter in the detector ensuring 
that only those electrons will be detected which have energy within 
a very narrow fixed range of energy defined by the filter. Doing so 
one gets NCE.co^ to be directly proportional to the Joint density of 
states , i . e. . , 


NCE,a):5 oc n^CE:»n^CE-ha)3 Ca. 5. lOI) 

provided we again neglect the energy dependence of the optical 
matrix element. Thus photo-emission experiment provides a very good 
study of the density of states or the band structure of the system. 

DISORDERED SYSTEMS : 

In chapter I we have already seen that an extra problem is 
associated with the study of disordered systems : that is 

configuration averaging. Therefore, to get an averaged optical 
conductivity, the right hand side of the equation C2. 4. 14D has to be 
averaged . This means, in general the product of square of the 
optical matrix element and the Joint density of states has to be 
averaged. This averaging can be done in two different cases : 
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CASE I : If the optical matrix element is assumed to be slowly 
varying function of energy and depends weakly on randomness* it may 
be separately averaged. Thus the vertex correction is negl ected * thi s 
factor is taken out of the summation in equation C2.4. 14D and 
averaged using the Virtual Crystal Approximation CPickett et.al, 
19S3D .Therefore, for semiconductors, we get 


<o<ooD> 


Srre^ 

m*^0 


L 

0) 




S <J,CE:>> dE 

tf 


where < > denotes the configuration averaging. This shows that, 

like ordered systems, we again have 

CO <oCcoD> oc S <J., CE:>> dE 

V f 

In this case photoemission experiment can again be interpreted in 
the same way as ordered systems. 

CASE II : In many instances is neither independent of energy 

nor randomness and its random variation may be correlated with the 
delta function part. We shall see that it is fairly strongly 
dependent on energy in dilute alloys, where the constituent 
potentials are very different, because of strong scattering by 
impurity sites. Hence the assumption that has been made in the case 


I is invalid. 
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In general for any three physical quantities A, B and C 
which have a relation of the form of A = B C , we have 

<A> = <BC> <B><C>. 

Therefore in such situation we have to take the average of the 
product of the square of optical matrix element and the Joint 
density of states in the right hand side of the equation C2.4.142> 
together. This means we have to include vertex correction also. 

<oC0)3>= - i- / r |g., |V,Cl-r,36CE -E. N C2. 5. 121 

Thus, in this case the energy distribution of 
photoelectrons NCE,a)D in the photoemission experiment does not come 
out to be directly proportional to the joint density of states, 
although very sharp features in the density of states will be 
reflected in NCE,wD, for example, neither weight nor the exact 
centroids of such structures will quantitatively match, due to the 
dependence of optical matrix element on energy and randomness. ; 

This is the case which has been considered in the present 
work and the averaging is done using Augmented Space Formalism. 





CHAPTER III 


AN AUGMENTED SPACE FORMULATION OF OPTICAL ■ CONDUCTIVITY FOR 

RANDOM SEMICONDUCTING ALLOYS 

3.1 INTRODUCTION : 

It was shown in Chapter I that the configuration averaging 
is central to the problem of the study of any disordered system. 
This is because the system has a statistical behaviour due to random 
potentials. It is also seen that it is the configuration average of 
physical observables which reflects the macroscopic measured 
properties of the disordered system. As we have already seen in the 
case of optical properties of such systems that, to study the 
averaged conductivity, we have to take the average of the product 
l%fyl -hcoD , where is the optical matrix element, i and 

f refers to initial and final states respectively and y refers to 
the direction of polarization of the electromagnetic wave. 

The Coherent Potential Approximation CCPAD introduced by 
Soven C19673 and later developed by Kirkpatrick et.al. C1970D , 
Stocks et.al. C19713 etc. provides a very good study of averaged 
density of states. Soven introduced the idea of a complex site 
energy Cor potential^ which characterizes the effect of the random 
environment within a meaningful picture. This potential C cal led the 
coherent potential^ was determined self-consistently from the 
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requirement, t-hat, a particular configuration immersed in the 

effective medium Cthe mean fields) yields no extra scattering on the 

average. Thus all scattering from single potentials, and 

uncorrelated multiple scattering from potentials are taken into 

account in CPA. But this approximation fails to provide the study 

z 

of the average of a quantity such as <5C -E^ -hcoD . Thus one 

needs to resort to a direct way of carrying out the averaging. 


The Augmented Space Formalism provides the best way for 
this propose.lt is a general method for the configuration averaging 
of any general function based on the direct way of averaging. Thus, 
it enables us to get an expression for the optical conductivity 
without making further approximation while doing the averaging i . e. , 
it includes the vertex correction. 


3.2 COHERENT POTENTIAL APPROXIMATION : 

The electronic density of states per atom for any system 
is directly related to the single particle Green function for that 
system as follows 

nCE> = - - ^ Im Tr G CE+iOD C3. 2. 15 

n N — 

where N is the total number of atoms in the system. The Green 
operator Cor resolvents G is given by the Hamiltonian H of the 
system : 

- HS”* 


GCzS 


Czl 


C3.2.25 



36 


where z is a complex variable. 


Therefore, for any disordered system the electronic 
density of states will be given by the configuration averaged single 
particle Green function. In case the disordered itself is 

’homogeneous’ - that is, the variation of the potential is similar 
throughout the system then the electronic density of states will 


be: 


nCED = ^ Im Tr <G Cr , E+iO:)> 

TT N — 


C3. 2. 3D 


where r is any site in the system. 

One of the main problems involved in the study of 
electrons in disordered systems is the study of configuration 
averaged single particle Green functions. One of the most powerful 
mean field or effective medium like approximation for the 
calculation of the averaged Green function is the self-consistent 
Coherent Potential Approximation CCPAD. The single site CPA had 
remained till late the main theoretical development in the study of 
electronic structure of random systems. It maintains the analytic 
features of the exact Green function and interpolates correctly 
between several individual limiting cases, e. g. , Virtual crystal 
approximation, atomic and dilute limits, in the case of binary 
al 1 oys . 


If we classify various literature on the single site CPA 
we find that these belong to three different C eventually equivalentD 


approaches: 
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Cl 5 Diagrammatic approaches using propagator or locator 

formalism introduced by Yonezawa and Matsubara C1966D , 
Yonezawa Cl 968:) and Leath Cl 968, 1970:). 

C2) The effective medium approach introduced by Soven C1967Z) 

C3) Multiple Scattering ideas first introduced by Anderson and 

M^'Millan C1967). 

All the different approaches yield the same final result 
within this approximation. This indicates that single site coherent 
potential is unique. Here only the second approach, i.e. the mean 
field approach will be discribed as it gives a very simple and a 
clear physical picture of the model and will form a basis of 
subsequent generalization of the CPA. 

The concept of a mean field or an effective medium in a 
disordered system was first clearly stated by Soven C1967). Seeds ; 
of the idea already existed in several problems on emulsions stated 
by Landau and Lifshitz. The idea of making it self-consistent 
apparently was Seven's. The trick is to replace the actual random ' 
potentials at different sites by an effective non-random potential I 
in such a way that the Green function corresponding to this i 

i 

constructed homogeneous potential gives the averaged Green function | 
of the actual random system. Thus if H denotes the Hamiltonian of 
the actual random system, the configuration averaged single particle 
density of states will be given by 



<nCED> 
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- - Jr Tr <G> 
rr N — 


C3. a. 4D 


where 


<G> 


<CzI 


HZ) ‘> 


C3. a. S3 


If H denotes the effective Hamiltonian of the effective 
— ef f 

medium then by the definition of effective medium we have 


<G> 


Czl 




C3. a. 63 


From the equations C3.a.53 and C3. a. 63 it is clear that H „ may not 

— O'!! 

necessarily be Hermitian. This means the constructed effective 
medium is defined by complex potentials sitting at various sites Cor 
the complex site energies3. Though the effective medium does not 
correspond to a real system, its Green function is related to the 
configurational averaged properties of the random system. 


Let us illustrate the procedure in a single band tight 
binding model. Our actual Hamiltonian in the tight-binding basis 


H 


i 


i-i 


^ r V T 


C3. 2. 73 


Here i is an index Indicating position or a site, say at r^. P. and 
Tj,j are projection and transfer operators respectively, [iXi | and 
|i><J[. £^J^s, the diagonal or ’site* energies and the overlap 

Integrals may be random. . Within a single site CPA it is not very 
profitable to include effects of off-diagonal disorder, l.e. , 


disorder in the V’ s. Therefore we set V. , = <V. ,> , where <V, .> is 

ij ij ij ij 
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some mean potential. Thus the equation C3. 2. 7D reduces to 


H = 



+ E 

ij 



<Vij> 


C3. 2. s:) 


Now we define our effective Hamiltonian In the same 
tight -binding basis as : 


H .. = E 2: CED P + E T. . <V, .> C3. 2. 93 

i o Tj -■‘J 

where site energies, are In general complex and energy dependent, 
but non -random and Independent of 1. 


yield 


The problem Is now to find out Z^CE3, which will then 


= <^z3> 


C3. 2. 103 


where z3 Is given by 




»eff" 


and 


<GCz3> = <CzI - H3 S 

— — — a.v 


C3. 2. 113 


C3.2.123 


Z^CE3 is determined self -consistently. In the single site CPA we do 

We embed an exact site energy at a site 1 


this as follows: 
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within the effective medium ,as shown in the figure C3. 2. in such 
a way that in doing so we do not produce any extra scattering on the 
average. This condition can be written mathematically as : 




^ff 


= <G> 


C3. 2. 133 


where G^^^ is the Green function for the medium in which all the 
sites are occupied by r^CE3 except the site i which is occupied by 
The corresponding Hamiltonian of this medium is : 


„<i> 


= 


Cs.-Z 3P. 

V o — t 


C3. 2. 143 


Therefore 


g‘^^ = Czl - 3"" 

— — — eff 

= fzl - H - C£.- Z 3 P.l 

j_ — — eff t o — Lj 

r Cc - Z 3P .,-1 

= Czl - H „3 ^ jl - — ^ S 
- -•ff [ zl- 

or 

= G „ fl - G ,c c«,- Z 3 P. 3 I C3. 2. 153 

— — •ff I — V o — V I 


If we expand the right hand side of the equation C3. 2. 153, we get 
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Figure 3. Z. 1 


Single potential embeded in an effective 
medi um 
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,CiD 




G 

-«ff 




C3. 3. leD 


= G „ 

—off 


— ®ff 


[c£ -E :»p ]e .Ji-G „[c/: -z dp ]] 

L V O — L — «ff »- t O — tJJ 


“1 


or 


= G 


Cs . 


df 


+ G 


-E DP. 
o — V 


1-G _CC^. “Z DP D 

— ^ff V O —i 


>ff 


C3. 2. 17D 


Define 


C^. ~Z DP. 

V o — t _ ^ 

1-G -Z DP. D ~ 

f f t O —L 


C3. 2. 18D 


Then^ 




off 


^ff 


t. G 


«-f f 


Theref ore 




= ^ff 


G „ <t.>G „ 

— ©f f “t. — ©f f 


C3. 3. 19D 


Condition C3.S. 13D restricts ourselves to : 


<t, > = O 
— i 


C3. 3. 30D 


i.e. , 


< 


Cs. -Z DP. 

t o — i 

l-<G><CcrZ DP.D 

— I o — i 



0 


C3. 3. 31D 
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This condition in the final form gives us an equation for in the 
single site CPA. 

th 

For random binary alloy, the i site either can be 
occupied by A atom or B. Therefore, for a binary alloy without 
short ranged order, equation C3. 2. 21D reduces to : 


«a-2: 

A o 


l-<G>Cs-l. 3 
A o 


+ ci-c:> 


£_-Z 
B o 


1-<G>C£;^-Z D 
B o 


= O 


C3. 2. 22D 


where, c is the concentration of A atoms. This is the CPA equation 

for a random binary alloy. To solve this equation self -consistently 

a suitable initial choice G Cz3) is made for <G> to start with. We 

— o 

then calculate Z from the equation C3.2.22Z> and then from Z we 

o o 

calculate a new ^CzD. We keep on doing this until Z^ converges. A 
careful choice of the initial G^ usually leads to a rapid 
convergence. Care must also be taken to preserve the herglotz 
analytic properties CMookerJee, 1973!) related to the physical 
necessity that the energy spectrum be real and the density of states 
be positive, this must be preserved at all levels of iteration. 

3.3 AUGMENTED SPACE FORMALISM : 

The Augmented Space Formalism CAS=*D is a technique of 
configuration averaging where configuration fluctuations are 
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systematically taken into account. In the effective medium 
approaches the quantum dynamical behaviour of the system is 
described by a Hamiltonian suitably set up, whereas the statistical 
behaviour is, as it were, imposed from the outside by allowing the 
potentials Involved in the Hamiltonian to vary randomly according to 
some prescription. The Hamiltonian, in itself, does not describe 
the full behaviour of the system and has to be augmented with the 
distribution of the set of random potentials. In the Augmented 
Space method, the Hamiltonian, expanded to include the description 
of possible configurations of the system .completely describes the 
disordered system. It is not necessary to put in the statistical 
description by hand. In this sense, this alternative description is 
by far the most satisfying one. Since the statistical description 
is incorporated within the Hamiltonian, the configuration averaging 
is not an extra process, as in the effective medium approaches. The 
formalism is formally exact but, for practical calculations, 
approximations preserving constraints of physical origin may be 
generated. 


The ASF was first introduced by Hooker Jee Cl 973a, 1973bl> 

and later on, Kaplan and Gray C1976, 1977D provided a detailed 

exposition. It has been successfully utilized for averaging 
quantities like one particle Green functions, particularly when one 
has to go beyond the single site approximations Ci.e. , beyond the 
CPAD and has to include clustering effects, large off-diagonal 
disorder, short ranged order, etc. 
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In a substi tutionally disordered system the lattice site; 
are randomnly occupied by the atoms of the type A, B, C, ... .etc. 
So the site energy s. is a random variable. For £ or e- depending 
upon whether the site is occupied by an A atom or a B aton 
respectively. The Hamiltonian for such binary alloy acts on the 
Hilbert space 9C spanned by a tight— binding basis. It involves randoir 
occupation variables <N^> defined as : 

= O ; if i site is occupied by an A atom 
=1 ; if i^^ site is occupied by a B atom 

A set of values taken up by the random ocupation variable <N^> is 
called a configuration. A certain probability density p(<N^>3 is 
associated with various values of the configuration. We assume 
that the various N^s are statistically independent, i . e. , there is 
no short ranged order. In this case we have : 

PC<N^>:) = C3. 3. 1!) 

where p^CN^!> is the probability density of the Individual variables 
. Kaplan and Gray C1977D have generalized this to the case where 
various N£s are statistically dependent within a Markovian short 
ranged order model . 

Haydock Cl 9725 first noted that p^CN^D has the samei 
properties of positive definiteness and integr ability as the density- 


4e 

of states have, i . e. , p^CN^i> 5: O 

/p^CN^DdN^ = 1 C3.3. 35 

Hooker Jee Cl 9735 used this equivalence properties of Pj^CNj^5 and 
density of states, nCE5 and suggested that similar to density of 
states, p^CN^5 can also be related to some operator in a 

hypothetical space corresponding to the random variable , and 

spanned by a basis set <|f^>>, through the relation : 

p, CN, 5 = -- Im <f^ |[CN+i65I-M^^^]~*lf^> C3.3. 35 

*^i i n o' ' o 

Thus the quantities lf^>» M^^^and p^CNj^5 are respectively analogous 
to the eigenvectors, the Hamiltonian and the density of states. 
Subsequently Hooker jee Cl 9755 described the disorder field in 
analogy with a fer ml on/boson field and was able to give a 
diagrammatic representation of it in a way similar to Feynman 
diagrams for fermions /bosons. The disorder field contains within 
itself a full description of fluctuations due to disorder. Now here 
the problem is to find out an operator H^*"^ in some basis < 
corresponding to a probability distribution p^CN^5 so that equation 
C3.3. 35 is satisfied. This is Just a ineverse problem of density of 
states where an operator C Hami 1 toni an5 H is given and the density of 
states is found from a similar relation: 

n, CE5 = - - ^ Im <i 1 Czl - H5“'^ ll> 

X fl N * — — * 


C3. 3. 45 
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where n^CED is local density of states. 

The problem can be handled in a way analogous to the 
continued fraction expansion of the Green function CHaydock, 1972D 

G^^CzD = <i |CzI - HD“'^(i> C3.3.5D 

In order to find a representation of \ given the probability 

density p^CRD, one has to try to write it in the form of a 
convergent continued fraction Cthe path shown in the figure C3. 3. ID 
by full linesD of the form as shown below : 


Pj^ C £+i = 


•9m 


l/'n 

2 

£+ir)-a - b 

' 1 1 

2 

c+in-a - b 
2 2 


C3. 3. 6D 


Such an expansion is convergent only if all the moments of p^C^rD 
are finite. 

Having written the probability density in the above form, 
the operator can then be represented in some basis < |f^’'^>> as 

a tri diagonal matrix : 
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a 

1 


b 

1 


O 

O 


b O O 

1 


a b O 

2 2 


b 


2 


a b O . 

3 3 


o 


b 

. 3 


a 

. 4 . 


b . 

. 4 


C3. 3. 7D 


Let us solve the problem for a random binary alloy We 

consider only the diagonal disorder and nearest neighbour overlap. 
Then the Hamiltonian H for such a system in tight -binding basis is 
given by : 


H = E 


P, 


E E 

■ . j € ft. 


T. . 

■^J 


C3. 3. 8D 


where 

= |i><i| is the projection operator 
and ^ij~ transfer operator 


defined in the Hilbert space 9C spanned by tight -binding basis < |i>> 
and ft^ are the nearest neighbours of i . Here ej^s form a set of random 
variables. Neglecting short ranged order, the probability density 
for the <N^> can be written as 


p^CN^J = c6CN^-lJ + Cl-cJ6CN^D 


C3.3. 9J 
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'n [n^+1t7-1 N^+iJ ’ ^ 


n 


,. cCl-cD 

N, +i7)-Cl-cD 


or , 


p CN :> = C-l/rrD .^ni 11m = 

7?— >0^ N, +i77-a - b* 

' i ' 1 1 

NiiT9-a^ 


C3. 3. lOD 


where a = c 

1 


a = Cl -c> and = cCl-cD. 
z 1 


Thus the above continued fraction expansion of p^CN^5 shows that 
corresponding operator has rank 2 and has a 2x2 tri diagonal 

matrix representation given by : 




c ^IccT-c5 

4cCl -c2> 1 -c 


C3. 3. ii:> 


in a basis 




spanning the configuration space. 


Let us now consider the configuration averaging. The 
average of any function fCN^D, function of a single variable ► can 


be written as : 
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f = 


X fCNOp, CN, ? dN, 
- CD i i i 


C3. 3. 121) 


or , 


where. 


f 


-Cl/rrD 



+ + 
f CN, +i0 D g CN, +i0 D dN, 
i ^oo i i 


C3. 3. 13D 


,ci:> 


g 


M 
oo 


= (CN^ 


+iO 51 - M 


Ci^'Y 


-1 


C3. 3. 145 


Let us further make the assumption that fCz5 as a function of a 

complex variable z has no singularities on that part of a real axis 

which forms a branch cut of the function g Cz5 , or in a 

*^00 

neighbourhood of the branch cut. In that case we may write : 

fCz5 g Cz5 dz C3. 3.155 

oo 

The contour is taken around the branch cut of g Cz5 along the real 
axis and not including any singularities of fCz5 Csee figure 
C3. 3. 255. 

f ' i 

Since is a self-adjoint operator on * therefore 

equation C3. 3. 155 can be written as 



= -Cl/^ni5 


§ 


fCz5 <f I 
o I 


00 

f Cz-h5“* dpCh5 

00 


|fo> dz 


f 


C3. 3. 165 




Figure 3.3. 1 Representation of three mathematical form 

of a positive, definite, integrable density 
function with finite moment - 



FI gure 3. 3. 2 


The contour in complex energy around 

the branch cut of the resolvent. 



where 


g 


-CiD 


M 
oo 


Cz]) 


00 



dpChP 

z-h 


C3. 3.1' 


and pChZ> is the spectral projection operator of M 


rCii) 


Thus 


f = <f^| J fCh:> dpChD |fS 

or , 

f = <f^ I |f^> C3. 3. 18; 

where the operator has the same functional form of as 

fCzD is of z. Thus it is found that the configuration average car 
be expressed as a representation of a suitably constructed operator 
related to the probability distribution. The entire process can be 
generalized for the average of a function of several randoir 
independent variables. For such statistically independent random 
variables N^s we have 

PC<N. > D = n pCN.D 

V *. ■ I 

I 

and so the entire configuration space is defined by 

$ = C3. 3. 19D 


Where i^^''^is the space associated to a single random variable and 
the full Green function is then given by 


G <M 


,CiJ 


>] 




C3. 3. aOD 



Then the average of the function 



given by : 


C3.3.ai5 

or » 

^ C3.3.E2: 

where 

1^0^ ^ nf C3. 3.23D 


The configuration space S is spanned by the basis set 



C3. 3. 245 


where <|f >> spans the space 4>^ . 

i 

As we know that each space <^*'’'’cfor random binary alloysD is spanned 
by two basis vectors. |f^> and |f^> .the space § will be spanned by 
2*^ basis vectors. FC<^>5 has the same operator function of<M^> 

as FC<N^>5 is a function of <N^>. Finally we define 


w = n ® nf 


C3. 3. 255 


This space is called the AvL^^ttvertted. space. For a binary alloy. 
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• N 

consist-lng of N atoms, the Augmented space is spanned by N x 2 

basis vectors. 

Now we are in a position to write the Hamiltonian of anj 
disorder system in Augmented Space. Let us considder a general 
situation where both the diagonal and off-diagonal disorder are 
present. The multi band Anderson Hamiltonian in the Tight— Binding 
basis in the space SIC , for binary alloy, is given by: 


H = 


E E 

i n 


'in —in 


E 


E 

n, m 


V T 

in,Jm -in,Jm 


C3. 3. 26D 


where 


s 


in 




Cl-Nj^D 


and 


V, , = N, N, + V®® CI-N.DCI-NP 

in,Jm nm i j nm i J 


+ <N, Cl-Np + N,C1-N, 5> 

nm i J J i 


C3. 3. 27D 


C3. 3. 28D 


i,J are the site indices and m,n are the band indices. 

Here the only statistically independent random variables 
are take the values either O or 1. 

probability density associated with each N^ is given by : 


Thus , the 
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PiCNj,!) *= c 6CN^-15 + Cl-c5 6CN^D 


C3. 3. 39D 


where c is 'the concentration of A atoms. Equation C3. 3.26!D can be 
rearranged as follows : 


where 


H = + E E <5£ N, P, + Z E v'*’ N.N. T, 

- -B n 1 -in nm i J -In.J 


m 


,< 2> 


+ E E V CN, +Np T 


i i*! j n , m 


nm i J — i n , J m 


C3. 3. 30D 


S£ 


n 


- A B., 

C£ - £ ^ 
n n 


v‘ 

nm 


yAA yBB 

nm nm 


£V^ 

nm 


nm 


yAB ^ 
nm nm 


and 


C3.3.31D 




EE^^P.^+EEEEvf® T 
J n —in ^ Pjjjj — 

in i n J m 


CiD 

The M for each value of N^s has representation 


M 


,C i5 


c ^JcCT-cT 

4cC1-c5 1-c 


C3. 3. 323 


in the basis lf'’>, |f''> spanning <p . The Augmented Space method 

* O ' 3t 

then yields the Hamiltonian H in the Augmented Space as : 



sc 


H = H^®I+5:£p, ®M 
— — B — n— in — 

1 .n 




+ F F T 

.*-‘. nm — in,Jm — 

m,n i , J 




+ r v‘ ^ ^ F T 

^ nm " -In.Jm 
n , m i i* j 


_ ^./l) „ _(i> . v,<j>^ X 

® CM ®I +M ®I D 


C3. 3. 33D 


^ <8> 

H e S«f ® fi® 4 > 


l‘''\ in the above equation indicate identity operators in 

(Jc ) 

all subspaces <p except those superscripted. 


3. 4 CPA THROUGH AUGMENTED SPACE FORMALISM USING GRAPHICAL METHOD 


Having described the Augmented Space formalism we are now 
in a position to describe CPA in the Augmented space. We shall 
describe it with the use of the graphical technique CHaydock 1972, 
ThesisD because it provides a clear visual illustration, 
particularly when we generalize to the cluster CPA’s. The standard 
algebraic formulation may often hide essentials of the approximation 
in the complexities of the algebraic equations. 

Let us first describe the graphical method. By a graph we 
shall mean a set of ’Vertices’ connected by ’links’. These vertices 



and links have one lo one correspondence with a countable basis od 

representation of an operator. If M is any operator represented 

the basis < ji>> then to each element in the basis < |i>> the 

corresponds a vertex v^ . Associated with each vertex v^ is t 

contribution [ l/<i|M|i> ] to the inverse of M. With each pair 

|i>, |J> is associated a link 1^^ and a contribution <i [M|J> to t 

inverse of M. The link 1^^ may be directional meaning may m 

be equal to Cshown in the figure C3. 4.1D. We shall define 

path in a graph of length n as a sequence of n+1 vertices connect* 

by n links. It is denoted by P Ci,J,k,...,nD where i , j ,k , . . . ,m ai 

n 

vertices. The contribution of a path of length n is given by 


kCP^Ci.J.k kD] 


n+ 1 


= n 

i = i 



n + l 

n 




M. . 

-^3 


C3. 4. 1 


Now we can associate the representation of the Gree 
operator Gin the basis < |i>> to the corresponding graph. The Gree 
operator or the resolvent defined as 

G = Cz^ - C3. 4. 2 

then the representation of G is given by 

00 

G, .CzZ) = E E kCP 3 C3. 4.3 

— d 1 " n 

^ r» = 0 F es 

n n 


where is the set of all the paths connecting the vertices v 
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and V 


J' 


It is clear that for a general infinite lattice the number 

of paths between two vertices is too large. Any Indiscriminate 

approximation to equation C3. 4. 3D may seriously violate the analytic 

properties of the resolvent. The enumeration or even statistical 

estimation of all paths is not a tractable problem. In order to 

evaluate the Green function we go through the Feenberg 

renormalization procedure. For this purpose we introduce the 

definition of a non-inters&c ting path P^. A path = 

<v 1 V 1 ...V > is defined to be Tion-intersect tnrf if none of the 

o oi 1 12 n 

internal vertices v . . . v are the same. If v = v then the path 

1 n-i on 

is called a closed, or polygonal non-intersecting path.. 

We now gather together all paths starting from v^» going 
through v^ to v^ either directly or via other vertices. This 
renormalizes the vertex v^ Csee figure C3. 4.2DD. The contribution 
at v^ are all paths going from v^ to v^ but not touching v^. This 
is Just Green function calculated on a graph with the 

vertex v^ missing. If we do this for every vertex, we obtain at the 
end only non-interacting paths but the contribution kCv^D is 
renormalized to CHaydock 1972, Mookerjee 1979D 

k’Cv D = 0 , 1 , 2 ,. .. ,n-0 C3. 4. 4D 

n nm 

The superscripts denote that the Green function is calculated on a 
graph in which the superscripted vertices are absent and 



fA* f 


1 



<iltn|j> 


Figure 3.4.1 


Graphical represent-ali on of vertices and 
links associated with an operator M 


•- 

0 



H 

0 



Figure 3. 4. 2 


Representation of renormalized Vertices 



G. c z:) 

— i-j 


C3. 4. 


oo 

= E E k *cp D 

n 

n=0 P€S' 
n 

where Includes all the non— intersecting paths of length n from 
the lattice. This expression is identical to Feenbui 
renormalized perturbation expansion. Statistical estimation of tl 
non— intersect! ng paths is a more tractable geometrical problem. 

Now we shall describe CPA through Augmented spac 
formalism. The Hamiltonian defined in the Hilbert space Sf, can fc 
written in the site representation |n> as 


H 


nm 


= £• 6 
n nm 


V 

— nm 


The Hamiltonian defined in the Augmented Space V' = ^ 
written in the site-configuration state representation |nf> 


C3. 4. 6 

can b 
as 


H , , = m]? 6 + V 6, 

— nf,mf ' — f f' riTO — nm ff' 


C3. 4. T. 


where, is the operator in the space § associated with thi 

probability distribution and V is the operator in For a bimoda' 

distribution of the ^ for each n is given by the equatioi 

C3. 3. 32:> and the graph is a single link chain like figure C3. 4. 3: 

where a = 1-c, a = c, b = JcCl -c5 . The configuration averagec 
1 2 

Green function for this model is then given by the resol veni 
corresponding to the Hamiltonian H, i.e. 
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<G CzD> = <nf |CzI-HD~* mf > 

iiin o * o 


C3. 4. 8:> 


here, the problem is reduced to that of an ordered system with a 
much amplified set of eigenfunctions. 

Until now all equations are exact. Approximations are 
made while employing the graphical technique. Haydock C Thesis, 
1972D showed that geometrical approximations on graphs involving 
delinking do not violate Herglotz analytic properties of the 
corre-sponding approximate Green functions. The configurationally 
averaged Green function <GCzD> can be determined by considering 
self -avoiding paths between |nf > and Imf > in the full Augmented 

IQ ' O 

space V'. Walks in the Augmented space imply that an electron at a 

site labelled by n and configuration state ~ caLn 

be induced by H either to make spatial hops to one of the neighbours 
of n with the matrix element V, keeping the configuration state 
same, or it can remain on the spatial site while the configuration 
state at the site n changes according to the configuration state 

at all other sites remaining same. 

Starting from the vertex |of> thefirst steps possible are 
to the near neighbours 1^ in the Hilbert space X' as shown by single 
lines in figure C3. 4. or to a different configuration state of the 
site ’o' labelled by, say, Cof^:> by a hop shown by double lines. 

If R Cz:> is the contribution of all self -avoiding 
o 

ion-inter sec ting paths from vertex O and back in the space Sf, the 
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Green function in an ordered system would be 


g CzD 
oo 


z-R CzD 
o 


C3. 4. 9D 


In three dimensions, R CzD will include the 

o 

contribution J] V^g^^CzD , where i are various neighbours of ’o’ and 

g^^^ corresponds to the Green function calculated from a sub-graph 
in which the vertex ’o’ is missing as well as the contributions of 
all closed self -avoiding loops in Sf. In the Augmented space the 
configuration averaged Green function can be written as 


<GCz3> 


1 

z-RCz5-TCzD 


C3. 4. lOD 


where RCzD is the contribution of all self-avoiding paths from the 
vertex ’of’ and back that are either Ca.2> entirely in the spatial 
part of the Augmented space or Cb:^ self -avoiding paths in the 
Augmented space which includes configuration state— hops but do not 
form closed loops. TCzD is the contribution of all the closed 
non-interacting loops from ’of’ in the Augmented space. Figure 
C3. 4.5;> shows a self-avoiding closed loops involving both spatial 
and configuration state-hops. 

There are numerous loops like this and it is almost 
impossible to account for all of them. Therefore the essential 
approximation is now made here to avoid this difficulty. Our 
approximation is to delink all those closed paths having both 
spatial and configuration state hops, i,e. 
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TCzD = O C3. 4. li:> 

Once we have made t-his approximation our approximate delinked graph 
consists of essentially that in the ordered system together with an 
extra configuration state hop with the off-diagonal element of 
in the tridiagonal representation Cb = -J cC 1 -cD in the case of binary 
alloys^ as the link contribution. If starting from a vertex ’nf’ 
we hop to a vetex *nf^* then the sub graph with the vertex 'nf* 
missing is an infinite graph which is exactly similar to the 
original one. A direct consequence of this property is : 


e 


cof:> 


of , of 


G . ^ = <G > 

of , of oo 


C3. 4. lay 


This will not hold in the absence of the type of delinking adopted 
above. 


This approximation gives exactly the ICPA result. We can 
show this as follows : 

Let G^^Czy be the Green function corresponding to the delinked graph 
in the Augmented space and let 



V, <5 , ^ ^ 
i J . i +1 


C3. 4.13D 


C 

be the tridiagonal representation of the operator M 
subspace of 5, then 


gDL 


CzD 


1 

z-R cz-w^G czyy-v^G czy 

O 11 11 


i n each 



6 ! 


= g Cz - W*G CzD3 

OO 1 1 


C3. 4. 14: 


where 


e^cz. 


1 

z-R Cz-W^G CzD>-wf^ G, CzZ) 
o 11 i+i i+i 


C3. 4. 153 


In the conventional ICPA formalism the self energy Z Cz3 is defined 

o 

by 




Cz3 


1 

z-Z Cz3-R Cz-Z 3 
o o o 


= g [z-Z Cz3: C3. 4. 163 

®oo o 

Z Cz3 = GCz3 C3. 4.173 

O 11 


W^= 6s VcCl -c3 for binary alloy 

Thus equations C3. 4.143 and C3.4.163 are equivalent. This 

equivalence has been discussed in detail by Bishop and Mookerjee 
Cl 9743. 


3.5 CCPA THROUGH AUGMENTED SPACE FORMALISM 


The generalization of CPA to CCPA involves a cluster of 
finite sites embeded in an effective medium, thus including the 


effects of correlated scattering from a cluster. 


This 





approximation, apart from retaining contribution of all paths 
described in CPA, also includes the self -avoiding, non-intersecting 
paths which involve hops that change the configurations over the 
n-sites of the cluster. Such polygonal self -avoid! ng paths take 
into account of the multiple scattering by the cluster. 

The vertices and links CbondsD have to be renormalized 
accordingly. To illustrate the procedure we shall, as an example, 
examine the case of a aCPA: a 2-sites cluster is immersed in a self 


consistent 

medi um. 

If 

we del i nk 

all 

pol ygonal 

paths 

linking 

configuration changes 

of 

more than 

two 

chosen sites, we 

have an 

octagonal 

decoration 

of 

sel f -avoi di ng 

pol ygonal 

paths 

in the 


Augmented space, graph shown in figure C3. 5. 1Z>. The eight vertices 
of the octagon correspond to the two sites and their four different 
configurations C2x2^:). This is the rank of the subspace in 
Augmented space spanned by the cluster and its configurations. The 
octagon belonging to a bond has no connection to that of another 
bond through links in the Augmented space, which effectively ensures 
that the correlated scattering from three or more sites belonging to 
different bonds is ignored. The final effect is that each bond is 
renormalized as shown in figure C3.S.3D. This gives rise to a 
matrix self -energy : 

C3.S.15 

The aim is to calculate Z self -consistently. It fully 
describes the self-consistent medium in the 2CPA. 
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The bond renormalization procedure follows in two steps : 

Step 1 : 


The decorating octagon is not an isolated octagon 
involving only the two spatial sites CO and ID. Let us 
take the vertex A ,£.e. , in the notation of Mookerjee 

C1973D. The bond AE Clf^ to of^D belonging to the octagon 
is only one of the z bonds emanating from A to the z 
nearest neighbours of 1. The other remaining Cz-ID bonds 
Clf^ to 2f^, 3f^,...D also hang on to the site A. These 

bonds themselves in turn have their own octagons decorating 
them Cfigure C3. 5.2DD. Thus the bond AE is itself immersed 
in a self-consistent medium. The same is true for the 
bonds FG and HK of the octagon. We may take this into 
account by saying that the medium renormalizes the bonds 
AE, FG and HK. Once the renormalization is accounted for, 
the octagon is effecively isolated. 

The renormalized vertices a CzD and links a C.-z') are 

o 1 

found as follows : We first divide the whole lattice into 

two sublattices, namely CID and C2D. CID is an 

unrenormalized bond AE. It has a Hamiltonian : 

= V C T T D C3. 5. 2D 

— 2 — AE 

r 

C2D is a lattice £ which is the original lattice mintis the 
link AE and in which all bonds and sites are renormalized 
by E. It has a Hamiltonian 
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G(H,.) 


V, = 2CV^ 4 C'V^ i Vpp 

V 2 = V'^ ^-i C (1 - C ) V ^ ■ ^-1 V F: a 


-«- Vj = 2 (1 - C ) V'^4 ( 1 - C) ' y^% V'e b 

->— X— X- \/rTr-T)[v'%-c) v"'] 
ycTi<)[v'"icv'"] 




c 


Fiqure 3. 5. 1 


C (1-C) V 


{]) 


/ci:r-c)e 
= eg4Cc 
= eg4 (1-C) c 

Octagonal decoration corresponding to ECPA 
in Augmer.ted Space 
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represents an octogonol — 

decoration 


Figure 3. 5. a Renormalization of site and bond of 

octagon in ASF due to Infinite medium 



Figure 3. 5. 3 


Final renormalized site and bond in aCPA 
for the octagon in Figure C3.5.1.D. 


the 



70 


= Z Cz3 Z P. + Z CzD r T. . C3. 5. 3D 

— o . i 

X J^A.E l , J ?<A,E 

represents the diagonal blocks in the partition 
of the Augmented space into a subspace ID spanned by the 
cluster and its configurations and its complement ED*. The 
off-diagonal blocks H and H = H are given 

by 

pCi.zD ^ cz:> r CT + T . + T. + T ,D C3. S. 4D 

— 1 . — lA — Al — tE -Ev 

U»«A.E 


The Green operator, , corresponding to the 

subspace CID is then given by projection of the full Green 
Operator, G, onto the subspace CID, i . e. , 

G^*^ = P G P 


|z^ - 


_ j|Ci.2Dt * C3. 5. 5D 


where 


G 


CzD 


zl - H 




-P 


C3. 5. 6D 


with P +P = I C3. 5. 7D 

— i — 2 — 

-Pv 

Here we have used the notation X ' to denote the inverse 
of the operator X on the subspace y. 


It is clear from the above equation that the effect of the 
rest of the lattice hanging on to the bond AE is to 
renormalize the Hamiltonian to + a where the 


self -energies 
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C3. 5. SZ) 


C3. 5. 9D 


C3. 5. io:> 


St-ep 2 : 


The next step is to work with the isolated renormalized 
octagon. This octagon is to be renormalized again because 
each of the vertices such as A,E etc. are connected to 
other vertices through either a hop that changes the site 
labelling or a hop that changes the configuration. For 
carrying out this final renormalization we again divide the 
octagon into two subspaces ; 

Cl!) a bond B COf, lf:> and 

C2D the rest of the octagon M 

The corresponding Hamiltonians are given by : 




C3. 5. IID 



E 

I , ji»=Of,4f 


H. . 

— tj 


C3. 5. 12:> 
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Finally we obtain the renormalized site self -energy Z^CzD 

and bond self -energy Z CzD as 

1 


Z CzD 

O 

= H® 

— OO 

E 

i^o 

E 

k^o 

rjBM -BM 

~Oj — jk — kO 

C3. 5. 17:> 

Z CzD 

1 

AfR 

= H + 

— Ol 

E 

E 

k^l 

H G Hf 

j — j k — k 1 

C3. 5. 18D 


Equations C3.5. 9D, C3. 5.10D, C3. 5. 17D and C3. 5. 18D together 
provide a sel -consistent set of equations for the 
calculations of the self -energy operator Z which defines 
the effective medium. 

The Green operator may now be calculated 



ESP. 
, o 




T. . 
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= E I + CZ /%D r V I 

O — 1 . . — i j ~ 


«■ - J 


Z I + CZ /VD h 
o — 1 — o 


C3. 5. 19D 


where 


h = T’ V T 
-o 


<^z3> = cz^ - 


C3. 5. 20D 


= fcz-Z DI - CZ /\Oh 1 

|_ o — 1 -oj 

(-1:] [h-v E7 ]i - 1.]" 


where 


= r-^f ] 


g C?D 


C3. 5. 213 


V and gCf3 = Cfl - h 3 

— — o 


C3. 5. 223 


Once we calculate gC?3 by methods available for 
translatlonally symmetric Hamiltonians, the equations 
C3. 5. 93 to C3. 5. 223 are the basic equations governing the 
2CPA. 


For the full cluster CPA, we first partition the Augmented 
space tfi into a ID which is spanned by all sites in the cluster and 
their configurations and its complement ID’ = v''^- Th® rank of ID is 
nx2^, where n is the size of the cluster and the system is a 

binary alloy. We now replace by Hamiltonian H* = ^03' 

effective H’^» which is translatlonally symmetric and has to be 
self -consistently calculated. 



Using the partition theorem 


Gjd * Czl - HD 

= - Ho> - H(D1D- 5^“’^ C3.B.e33 

' - H|D - £ 


Here 


!1id y. ^a> 


^ln>D* ^ ii5[d 


-ID _ ^ o -V 


^«eff^ = 


ib-ID 

IDID 


C3. 5. a4D 


We now partition the space ID on which H = +o; acts into one c 

which is spanned by CnfD and its complement c'. Again using the 
partition theorem 


G 


C 


Czl 





= zl - H. 




g'=’ch^„3 


««■ 


c 


C3. 5. 25D 


But the Augmented Space theorem states 







C3. S. 26; 



Equations C3. 5. 24D and C3. 5. 26D are the seif -consistent equations 

<e> = CzA - h A = r A C3. S. 273 

— — O — — — 

where A = h H C3. S. 283 

— — o — ef r 


may be easily calculated, even if A is not diagonal by using the 
modified version of the Recursion method with A taking the place of 
the overlap matrix. 


3. 6 OPTICAL CONDUCTIVITY FOR RANDOM SEMI COHDUCTI NG ALLOYS IN 

AUGMENTED SPACE 

We have already discussed the general formultation for 
optical conductivity for any crystalline solid in § 2. 4. In this 

section we shall now use the basic formula for optical conductivity 
and generalize it for random semiconducting alloys within the 
Augmented space formalism. 

The basic formula for optical conductivity is given by : 
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oCftD 




1 

<A> 


E 

iTr 


l<V'f |e^.Pki>|*fiCl-f^3 


6CE^-E^-hco5 


C3. 6. ID 


where 


m 

0 
to 

1 
f 

r 

p 


effective mass of electron 

volume of the specimen 

frequency of incident photon 

refers to the initial state C valence bandD 

refers to the final state C conduct! on bandD 

refers to direction of polarization of the 

incoming photon 

momentum of the electron 

occupation probability of initial or final 
states 


For semiconductors, at T = 0°K, we have 


f, = 1 and f^ = O 
i f 


C3. 6. aD 


Thus, the equation C3. 6. ID reduces to 


oCwD = 


E |<V-f|e p|v-,>i' 6CE^-E^ 

m»t ID ify 


-hwD 


C3. 6. 3D 


As before, we use the property of Dirac 6-function and find 


otcoD 


2ne* 1 


z ^ <>> 
m* D i f r 


E MVfle ■p|»'i>r J '5CE-E^3 <5CE-E^5 dE 


2j E 

S /“ l<*'rlVPlVl^ E*:e-Ej 3 aCE-E,3 dE 

®" “ O C3.e.« 
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Where the upper limit E^^is determined both by the band limits and 
h(o. But, by definition : 

40 * E D 6CE-E D 6CE-hto-E.D = n.CED n, CE-hcoD C 3. 6.51) 

1. JT ^ 1 I X 

where n^CED : conduction band density of states per atom or 

bond 

C E-hci>D : val ence band density of state per atom or bond 

The factor 4 appears in the equation C3.6. 5:> due to two possible 
spin orientations of the electron in the initial and final 
states. Therefore the equation C3. 6. 4> reduces to 


p 2 . O /s _ 

- E r •p1w.>|* — — n.CED n, CE-h«:> dE 

„2- CO ^ j i I ^ ^ ' .^-2 f i 

m* O y ^ 40 

o 

E 

& rs ^ 

^ ^ |e •p|V'^> n^CE> n^CE-hw5 dE C3.6.6D 

S m ^ o ^ 

We shall be dealing with disordered systems and shall 
interested in finding the configuration averaged o<<oD by averaging 
the product of the square of optical matrix element and the valence 
band density of states. Although, in general, the product of the 
square of optical matrix element and Joint density of states has to 
be averaged to include the full vertex correction, but subsequently 
we shall factor the initial and final densities of states. We shall 
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take the conduction band density of states to be free electron like 
in a Virtual Crystal sense. Since the conduction electrons are 
reasonably extended, such an averaging may not be too bad. The full 
calculation should include the Joint two particle propagator for 
initial and final states. The Augmented Space formalism is capable 
of handling such problems CThakur et . al . 1987:5 . We shall leave this 
development for later work. As the free— electron like density of 
states is randomness independent. 

We have, 

<oCc«>D> 

av 

C3. 6. 7:5 


~2 § ^ J { |<Wf.|e .plv/^>i^n^CE-hcoD \ dE 

2m* ^ O 


We shall expand the initial states in terms of the set of 
Tight-binding basis <[n>>. For a semiconductor the Ti ght -bl ndi ng 
basis <jn>> consiting of say, bond orbitals Cto be discussed in 
§4. 2Z> , do span the valence and subspace. Therefore equation C3. 6. 7Z) 
gives the averaged optical conductivity <o<i(jD> to be 


e^ 0 
2m 


S-m. 


E E E 

m n 


6 

mn 


J n^CED^<¥/^ jp^lmXml^p":) |n><n | p^ | dE 

O 


C3. 6. 8:) 


Here we have replaced the valence band density of states n^CE-hcoD by 
the corresponding Green function. 
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HiCEi - = - ^ Im E E :> |n> 

m n 

where 7 ?" = E-hw-O''’ 



We may expand the final state in terms of plane waves 


V 


f 


^iCk+G3 

G 


V' ■ 

G,k 


C3. 6. 9:5 


Replacement of the full sum over all reciprocal lattice vectors by 
only G = O leads to the free electron case. The optical matrix 
element factor in equation C3.©. 85 can now be evaluated as follows : 

<n\p^\^^> = <nle^.C-ih^ |V'^> 

= -ih<n|e^. C3. 6. 105 


In the r representation, <]n>> are given by the centired orbitals 

like <d)Cr - F 5>. Therefore, 

^ m 


<nlf> 



_ i 

r 5 e 
n 



C3. 6. 115 


Therefore equation C3. 6. 105 becomes : 



ih 


jAi 


/V i k^.r 

r 5 e . 7 e dr 

n y 
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ti 

w 


✓v 



r > 
n 


d r 


h 

w 




C3. 6. laD 


1C ^ 

where <p Ck^D is the Fourier transform of <^rD. But for random 
alloys ,<^r3 ’s are random. For a binary alloy of constituents A and 
B with concentration c and Cl-cD respecti vi ty, 0Cr — r^D is given by 

0CF-r ■>-<(>. Cr-F D N + <j&^Cr-r 5 Cl-N !) 
n n n ^B n n 

= <^^Cr-r > + 6<^r-r 5 N C3.6. 13D 

^B n ^ n n 


where N is defined as : 
n 

*1 ifn site is occupied by A atom 

= 0 if n*"^ site is occupied by B atom. 

0^CF5 and <^gCFD are respectively centred orbitls associated with A 
and B atoms and 64> ~ 4>^ ~ 

The probability density is given by 

pCN > = c 6 CN -ID + Cl-cD 6CN D C3.6.14D 

n n n 


Therefore, 



C3. 6. IE 


* _ ik 
4> Ck^D e 


f • 



[<^gCk^D + 6Ak^DNj 



r 

n 


Hence for random semiconducting binary alloys, the equation C3. 
becomes : 






C3. 


Si mi 1 ar 1 y ^ I ■ P j ni> is gi ven by 


<V'f |ej^- Pj»n> 






C3. 


Equation C3.6. 8D now becomes 


<„CC0» ^ S "fCEJ / [- 

* y rn n ' 


am 


o 


. /s -ik r 

hr- f m 

e. k^e 

^5 


<fl!.gCkj.D + 6.#Ckj.5N^>1 <m|K7)”3 |n> 1^ 

-i ^ jj 




e . k^ e 

r f 


ik„. r 
f m 


<<P^Ckp + 6<p'^Ckp 


v]> 


dE 


O 


2 CO 


2m 



y m n 


mn 


E 

J 

o 




6 . 12 


6. 16; 


6. 17; 


' D 
. m 
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+ a^k^NjG^CT,-3 [,^*Ck^5*6Al=f3Nj^ dE 


^ I E E E <5^ J °nfCED < 

c2m y m n ^ ' 


ik^. Cr -r D 
,2 f n m 

k^ e 

tr 


\<p^k^:) ^ G Crr :) + Ck.D 64>Ck^':) N G Ct) D 
Bf mn' f mmn' 


<p^Ck.:>6<p CkpG Crj DN + j6<?l>Ck.D *N G Ct? :)N ) dE 

^B f ^ f mn ' n ' ^ f m mn ' n | / 


C3. 6. 183 


wher e k , 


e . k^ 

r f 


Equation C3. 6.183 clearly shows that to average <yCoD we must average 
the following quantities : 


C13 

1 vv 



C23 

**Ck^>6^k^3N^G^C„-3 



C33 



C3. 6. 193 

C43 





ik .Cr -r 3 

The factors k^ and e are randomness independent. 

tr 


The averaging will be done using the Augmented Space formalism. Let 
us first write the averages of above quantities in the Augmented 



Space with site configuration state representation : 




csy 


C4D 


' L J 

= <mf I ||nf> 

' L Jav. 

= <^gCF^3<5<j!Ck^D <mf|^ ® C3. 6. 20D 

= ip^Ck <Jiif |G ® |^jnf> 

/ < |6^kp I^N G ^C7 ?"dN^\ = |60Ck-D|*rNG Ctj'dn] 

\ ' f' mmn' n/ f'|_mmn' ’^Jav 

|6<;(!Ck^3 |* <mf (^ ® G ® ^|nf> 


We shall now evaluate the above quantities using the representation 
of ^ in configuration space <p^, given by the relation C3.3.32D. 

Therefore : 

<mf (]^ ® G|nf> = 

+ T 11® G|mf> . 

mm! — * 

Y y J 
^ O 

Here we have used the notation for 2 ind y™ for |f^ in the 

configuration space <p^. 


+ Cl-cD 


cCl-c> 


m m 

r Y 

a * i. 
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p = |f"V<f'"| 

m ' o o' 


and 




<mf 1^ ® e]nf> 
<mf |G <» 1;^ jnf> 


= c<mf|G|nf> + 4 cC1-cD <inf^jG|nf> 

|G <» fcP + Cl-c^P 
'—In n 

L V y 

O 

+ 4cci-c:) Ft + t 1] |nf> 

Inn n n I I • 

*- V V V V “* -■ 


<mf 


y y 

^ O* ± 


T T 
* ±* a 


= c<mf G nr> + icCi-cD <mf |G nf > 
— I— n 


and 


+ Cl-cD P 
m m 

y y 

* o * ± 


<mf Im"*® G ® M"|nf> = <mf‘| fc P 

^ 7 

+ 4cci-cD Ft +t 11g®FcP 

^ I m m m ml — I n 

L. y y y y J J JL y 

* - T ^ ^ ] 1 |nf> 

L V V V V J J 


+ Cl"”cD P 


n 


y y 

^ & t 


r r 

1 o 


l^c <mf| + ^cCl-cD < i j^c|nf> + 4cCl-cD inf^>j 


c <mf |G nf> + c4cCl-c<inf Gtnf> 
* — m — 


+ c4cCl-c3<inf iGlnf > + cCl-c3 <mf Glnf >. 
^ 1 — 1 n m — I n 


Thus the expressions C 3. 6.205 now become 


C15 

C 25 

C 35 


4>QCkp |*<mf |G|nf> 
gC kp 6ipC kp j^c< mf | 

^gCk^5<5<^’^Ck^5 |^c<mf |G|nf> + JcCl -c5<inf iGjnf j 


^gCk^56^k^5 j^c<mf |G|nf> + ^cCl-c5<mf^|G|nf>'] 


C3. 7. 215 
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C4D 


I* |Gjnf> +c-^cCl-c3 <mf | Gjnf > 


+ cjcCl-c5 <mf iGlnf > + cCl-cZXmf iGlnf'^- 

* — ' n )xi * — ‘ nj 


Therefore^ after taking together the above expressionsC3. 7. 21D , we 

get 

®j<mr |Gjnf> 




B 


t- 4 cC1-cD 0g00 + c|cCl -c:> j 60 1 ^ j < mf ^ j G j |nf> 

h |V cCl-cD 0g60** + cV cCl -cJ> j 60 1 * j < mf | G j nf + 

cCl-cT) |60|^ <mf^[Gjnf^> 

<mf |G|nf> j^0g0g + c0g0^ - c0*0g + c0g0* - c0g0. 


+ <mf iG nf> 
m' — 


c 0A^A ~ *^A^B “ ■= ^a'^B 


4cC 1 —cD ^cOr*~cT <P^<p^ + c^'c'C '£™ q3''' 4^ .<P j. 

o A rs r> A A 


+ c4cC1-c 5 0g0g- c^cCl-cD 0^0* - c4cCi-cD 0*0gj 
+ <mf |G[nf^> |^^cCl-c5 0^0* - TcHT c^|— 

+c^cCl-c30g0g - c^cCl-cZ) 0^0g - c4cCl-40^0gj 
jGjnir^> |cCl— cD6060 J 


<mf 


<mf iGlnf> [^0g0g 


Cl-8c+c^5 + 0^0* + 0^0gCc-cb+ 0*0gCc-c*:> 


] 


<mfj^|G|nf> j^0^0g f/cCl -c3 - c4 cCl -cD ) -0g0g Q cCl -c5 
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Ci|cCl -c^) + c4 cC 1 -cD ~ c4 cC 1 -cD J ■*" 
■'■ I ® I l^c^cC 1 -cD 4>p4^ ”0g<^g CTcO~--c 5’ - c4cCl -cd} 

- c4cCl-cD<^^<^g + Q cC 1 -cD - c4cC1-c33 j 

+ < mf ^ j e I nf FcC 1 -cD 


- 


<mrjGjnf> 1 0g<^gCl “cZ) ^ ^^^gCCl —cZ) + 


[^4cC1-cD - <^g<^gCl-c5 + 


0>bCcD 


]] 


+ < mf I G I nf ^> 1^4 cCl -cD 
+ <mf^|G|nfj^ j^cCl-c56<^<?!>*J 


<mf iGlnf> jci-c3<^*<:Cl-cD<^g+C(^^> + c <^*<c <^^+Cl -c><^g> j 

+<mf^jG|nf> CcCl-cDD^''*[^Cl-cD0g ■*■ ‘^A^^a'^B^ 

+< mf I G j nf ^> c cci -cD [ci -cD 0gC <i>^-<P^ + c<p^c ^*-<^g:) ] 

+< mf I G I nf > cCl -c]> 6<p6(^ 

m ' — ’ n 

<mf |G|nf> j<c<^^+Cl-cD<^g><c<^*+Cl-c5(^g>j 
+< mf ^ I G I nf > C cCl -cD D |c <c<^^+C 1 -c5 ^g> j 
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+ <mf le|nf^>CcCl-c5D‘^* -cD^g> j 

+ <inf^|G|nf^>cCl-cDC6^<;«>*] . 

But c<p. + Cl-cZ)<;frj^ = ^ 

A B 

c<^* + Cl-c5<;frg = ^ 

~ = ^4> C3.6. aa? 

<?!>^ - <!feg = 6<p and 

<mflG|nf> = EG ] 

' — ' mn av. 

Thus the sum of the four quantities C3. 6.19;> becomes : 

S = EG ] E0. ^3 + <mf I G I nf > C cCl -cD 5 L6<p^l 

mn av. ^ ^ m' — ' ^ 

+ <mf jGjnfj^> CcCl-cDD*'^^E6<^*^] 

+ <mf iGlnf >cCl-c:>E606<^*3 C3. 6.335 

m'— ' n ^ ^ 

where the sum of the four quantities in equation C3. 6. 195 is denoted 
by S. 


Since the Augmented space is homogeneous 
using the symmetry of the Augmented space, we have : 


therefore by 



C3. 6. 345 


Then, 
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^ <Jnf',j,|e|nf>CcCl-cDD*^*C 6 ^+ 6^*'^3 

But., 6<p4> + 64> tp - ‘2, Re C^* Sp"^ , therefore. 


S = Ce^^] ^ \p\^ + <mf jG|nf>CcCl-cD5^ 8 R&Cp tp:> 

mn av. * » jn» — * t- 

+ <mf^|Gjnf^>cCl-cO j6<^|^ C 3. 6. 855 


Let us now evaluate <mf jGlnf> and <mf |G|nf > using 

m ' — ' m ' — ‘ n 

graphs. We shall use different notation now that is 



®mf , nf 
m 


<mf G nf > = 
m — 


®mf , nf 
m n 


In graph G^ nf given by 
m’ 


A# 

G 


inf , nf 
m 


2 


kCP5 


C3. 6. 865 


PeS 


mf , nf 
m 


where S ^ ^ denotes all the self -avoiding paths between vertices 

l!!lX f> iTliX 

m 

mf and nf. But can be written as Csee figure 3.6.13 

m mf ^ , nf 

m 


P - p 

^mf ,nf ’’ *^jnf.nf *^mf .mf 
m mm 
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where 4cCl cD is the contribution associated with the link 1 ^ ^ 

mi , mr 

m 

and 


~ and Cg are the site energies 

corresponding to A and B atoms respectively. 


Ther ef ore ®mf , ^ can now be written as 

m’ 


G = G k IP' 

mf „ . nf mf , nf '• mf , mf 

m mm 




But 


(pmf 


.mf 3 

m m 


m m 


Hence 


5:mf 


.nf = ®mf,r,f 

in mm 


Fg 1 ^ 

I mn I ^ 

k Jav. 


cCl-c!) G 




mf , mf 
m m 


C3. 6. 27:) 


Similarly, we evaluate G^ Csee figure C3. 6.255, i . e. 


m n 


J kCP__ .. 5 


mf . nf 
m n 


mf , nf 
m n 


P = p”*^ P P^^ 

mf , nf mf , mf mf , nf nf , nf 

m n mm n n 


So, 


G ^ ^ = JcCl-c5 6s G^ ^ O ^ G*^!! „ .JcCl-c56£ 

mf ,nf ^ mf^.mf^ mf,nf nf_,nf_ ^ 

m n mm n n 


C3. 6. 285 
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Figure 3.6.1 Graphical representation of the path 

P , Iru Augmented Space 
inf 

m 


tn4 

/#/ #- 



vt’Ci-C) 


m 



<a~aphical representation 

p in Augmented Space 

mf ,nf 
m n 


igure 3. 6. S 


of the path 
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= EG 3 

mn av. 


C<5^-^cCl 


Putting these back into equation C3.6. 251), we get : 


S = EG ] 1^1^ ^ ^ i| cC 1 

mn av. * ^ * mn av. 


—cy 


2Re C^6<p:> + [G _3 2 

CSeicCT^^^- 


cCl-c3 \ 64>\' 


IG 3 
mn av 


1^1* + 2Re 


Sc 


l<S<>i^] 


CG 3 
mn av. 


[l^ + 


Calling p + ^ Sp by finally get 


S = [G 3 , 

mn av. "^eff 


C3. 6. 30D 


Therefore the configuration averaged optical conductivity for random 
binary semiconducting alloys is now given by : 


* 2 E ik..Cr^-r > 

®-iI- - S-tw EEE <5 f nCE D e 

2m o 

CG 3 |*clE 

mn av. '^eff. ' 


< ad <)D> «= 
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e*h*n 1 „ r ° 2 . 

^5 L S n^CED 

o 


. ik-.Cr -r 3 , 

^ Z E 6 e ^ ™ [G ] dE 

” S n mn av.J 


2 j, 2 ^ r 

ne fi 1 ,2 


- - .<Pflt r CG 3 dE 

n " mm av. 

L Tn J 


Therefore 


<cxCw>> 


2^2 
Tie h 


I CE kj 3 CED<n CE-h«D>dE 

J y * 


C3. 6. SID 


Let us now evaluate S}c. , where k_ = k^.e . Choose the direction 

^ fr fr f r 

f 

of along the z-direction Csee figure C3. 6. 3;5D, e^ being the 

direction of incoming photon. Therefore, k^^ = k cos & or, k*^ = 


I 2 Z- 

k cos 6. 


For an isotropic wave we have E “ where, is the solid 


angle subtended at the origin. So the sum = J” dCcos 0D d*;^. 


Therefore 
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1 zn 


E k 


fr 


■J I*; 


COS B dCcos ©D 


“1 o 


% 

I 


errk^ f cos*© d cos © 


4 rr , 2 
-3 


Putting back this value in equation C 3 . 6 . 311 ) we get 


<oCo) 5 > = 


g^»t 2 w J 3 i*^eff I r^f^^EiXn^CE h 6 >D> dE 


— 2 2^2 
err e 


3iti 




kJ|<A^^^|*n^CE:)<n^CE-hcoD>dE 


. 2 , 2 

k 

But is related to E through ECk^D = — — — , therefore k^ = . 

^ Em* ^ h* 

Thus finally we have the expression for average optical conductivity 
for random binary semiconducting alloys as : 


<< 7 Cw 5 > 


.22 
4 tt e 

3 m 


o 

S J E^eff 


O 


.22 
An e 

3 m* 


o 

- I lg n^CE 5 <n.CE-h<o:>> dE 

(o J > ’’•f f ' f 1 


C 3 . 6 . 32 D 


0 
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where the square of the optical matrix element g CED is given by 


3. 7 APPLICATION TO A MODEL SYSTEM : 


In this section we shall apply the formulation of optical 
properties of random semiconducting alloys, described in § 4.6, to a 
model system. For our model we shall choose the following 
parameters : 



= -5.0 ; 



= -9.5 ; 


V 

= -1.0 

and 

c 

= 0. 5 and 

0. 2 


where and are the site energies for A and B atoms 

respectively, V is the overlap integral Cno off-diagonal disorder!) 
and c is the concentration of A species. 

The Hamiltonian for the system, in Tight -binding basis, is 
given by the equation C3. 2. 8>, where or and * V- 

The valence band density of states for the system is calculated 
using ICPA and 2CPA. For our ICPA calculation we shall choose 
rounded density of states Cshown in figure 3.7.13: 




9 ' 


r 

n CEO = i 
o i 

K. 


-| Cl-E'i*-"* 
0 


lfjE| <1 
if|E| >1 


C3. 7. 1 


as the initial choice for <G> appeared in the ICPA equatioi 
cs.a.as:). We then calculate self-energy E^CED using the ICPy 
equation C3.2. 22D and then new <G> self consistently. The valence 
band density of states <n^CE^D> for c = 0. 5 and c = 0. 2 are shown ir 
figuresCS. 7. 2D and C3. 7.3D respectively. 

The 2CPA calculation of <nj^CE^D> is done by choosing the 

initial choice for Z as calculated within ICPA and Z equal tc 

o 1 ^ 

<V. .>. We then calculate new Z , Z and <G> self -consistently usinc 
x j o 1 ^ _ 

equations C3. 5.9D, C3.S. lOD, C3. S, 17D and C3.5. 18D. The averaged 

density of states <n^CE^D> as calculated within 2CPA for c = 0. 5 and 
c = 0. 2 are shown in figures C3. 7. 2D and C3.7. 3D respectively. 

There are some structures seen in the impurity band in the 
case of 2CPA calculations of density of states. These two peaked 
structures are due to bonding-anti bonding states of the AB cluster. 
These structures were also seen and extensively discussed in the 
earlier works for both model calculations CKumar et.al. 1982D and 
for Cu-Ni alloys CThakur et. al . 1987D. 

We shall now calculate the optical conductivity, oCojD, and 
other related response functions. The optical conductivity is given 
by the equation C3. 6. 32D 

The conduction band density of states, n^CED, Cfor free 
electron likeD is given by 



Figure 3.7.3 
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n CED = rr^- 

C2n:)^h 


If we choose units such that h = e = m* = 1 , we have 


<o<.oD> = 


4n^ if 2 

"3- j;; ^l^eff I — ^ <nj^CE-h<o3>dE 
J *n 


03/2 if , 

— 5 ^I'^effl JE <n CE-<^>dE 


C3. 7. 2D 


From the conservation of energy we have 


E^ + CO = E. 

Since <n^CE-coD> is non-zero between the band edges E^’^'and 

and zero outside the edges Csee figure C3. 7. IDD, the integrand in 

the equation C3. 7.2D is always zero outside the range between e^^' 

TY ^ A "V* 

and E ' . If we change the variable from E to E. , we get : 
o 1 ” 


< <yC coD > = 


2®^* 1 
3 CO 




+coD \<p 


<n^CE^D>dE^ 


C3. 7. 3D 


Let us now evaluate given by 


E CE D 

^eff^^^ = ^ED + ^ <5<^ 

^ = c ^^CED + Cl-cD ^gCED 


C3. 7. 4D 
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~ ~ and 


A - ^B- 


We shall choose <fr^CrZ> and ^gCrO t-o be t,he wave functions C t=03 for 
a square well potential of width 'a' and depth CSchiff 1968D. 




\P 


A 

Stt 


-ft^r 


and 


= J 


2fr 


-V 


where I is “the angular momen'Lum quant. urn nuinber » 


2m I E . I 

and ^g = CajEglD^- 

^A are the eigenvalues of the respective atomic Hamiltonians. 

If we assume that the atomic potentials are not much afffected 

o 

within the Winger -Seitz cell when the solid is formed, within the 
first approximation E. and E^. may be taken as site energies e, and 

A ±5 A 

respectively. Thereore 

jy 

and ^g = 

The Fourier transform of <^^Cr3 and <^gCrD are respectively : 

JEftl _ 

^ Ck3 = and Ck3 = — ^ C3.7.53 

Now using the equation C3. 7. 32) together with the equations C3.7. 4D, 
C3.7. 53 we calculate the valence band density of states <n^CE^3> 
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and the self energy E^CE^>. We then calculate the optical 
conductivity for the model system. Then using the various relations 
among the response f unctions > descri bed in chapter II, we calculate 
the other related response functions. 

Figures C3.7. 40 and C3. 7. 50 show the optical conductivity 
as a function of frequency of the incoming photon for the model case 
with concentration c=0. 5 and calculated within the ICPA and the 
ECPA respectively. Figures C3. 7. 60 and C3. 7. 70 show the optical 
conductivity for the model case with c=0. 2 and calculated within the 
ICPA and the 2CPA respectively. Each of these figures clearly show 
the effects of disorder which arise mainly due to vertex correction 
associated with averaging of the kind of < |*. n^CE0> . The effect 

is seen in lower part of the frequency only which is expected 
because it is in this range of energy that [g.^^CEO j has 
structures. In the higher frequency range the factor 1/to of the 
optical conductivity dominates over others making it decrease with 
frequency and washing out any disorder effect. 

There are small double humps seen in the conductivity 
curve. These humps are the reflection of the double peak DOS 
structure and the structures in the square of the optical matrix 
element these humps are much broadened because of 

convulatlon of the valence band density of states with the 
featureless conduction band DOS and relatively featureless optical 
matrix element. The effects of disorder in the case of c=0. 5 is 
more prominent than the c=0. 2 case. This is expected because in the 
earlier case the "degree of disorder is stronger than the latter. 
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Figures C3. 7. 8:) to C3. 7. 155 show the other related 
response functions like real part of the dielectric function 
C related to permitivity5 , imaginary part of the dielectric function 
, refractive index and the absorption coefficient. Each of these 
figures shows the effects of disorder. The double hump structures 
are also reflected in the real part of complex dielectric function 
and the refractive index function. 

These figures show that the ICPA and the 2CPA results of 
optical conductivity and other related response functions do not 
differ much. They are almost the same. The structures which arise in 
the SCPA EXDS are not reflected in the conductivity because it is 
related to the integrated CconvulatedD density of states. 

Figures C3.7. 165 to C3. 7.235 show the square of optical 
matrix element |g I* as a function of energy for different values 

I ' 

of the concentration ,c, and the frequency of the incoming photon ^ 
for the ICPA and the 2CPA calculations. Each of these figures shows 
that the depends on both , the energy and the disorder. The 

|q I* calculated within VGA has a weak dependence on energy . 

t s'tfy I 

However , disorder effects lead to a stronger energy dependence of 
the optical matrix element . These figures show that the effects of 
disorder not only change the magnitude of also shows 

structures . These structures are more prominemt in the case of the 
2CPA and c=0, S . These are expected because in the case of CCPA the 
effects of cluster also come and for c=0. 5 the disorder effect is 
more . The structures in arise due to the self energy , as 

the correction term is given by CZ^/<5£5<5^ - It is also clear that 

£ A ^ - ■* - 
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because the self energy has sharp peak near the impurity band peak 
,t.G. at E=-5. 0 . This sharp structure. is not reflected in the cas« 
when o) -4. 6S because the transition takes place only from that 

portion of the valence band for which E =-4. 62 to O. O , whereas ir 

i 

the case of co =9. 24 , the transition takes place also from t,he peak 
of the band. 

Figures C3. 7. 240 to C3.7.310 show the product of valence 
band DOS and the square of the optical matrix element, 

^ function of energy, E, for different values of c 
and (0 with the ICPA and the 2CPA calculations . It is worthwhile to 
notice in these figures that , though the DOS for c=0. 5 and c=0. 2 
are very different to each other , but when the DOS is multiplied by 
the square of the optical matrix element , it so happens that the 
product <n^CE^0> |g*j|^CE0 look to be similar except slightly 
different in magnitudes . Since in an energy resolved photo ejnission 
experiment one essentially studies <n.CE.5> jg I as a function of 
energy , therefore the results of such experiments do not reflect 
the jlolTit' DOS but very much modified by the square of the 

optical matrix element which arises due to energy and disorder 

dependence of (g CED | . Thus in order to get .correct band 

t fy 

structure from such experiments , one must factorize the square of 
the optical matrix element factor . 

Though , as it has been seen above, the 2CPA calculation 
does not make much difference from the ICPA calculation in the study 
of integrated response functions , but the results of the square of 
the optical matrix element show that the 2CPA result differs 
significantly from ICPA result . Thus in the band structure 
calculations one must take into account the cluster effect. 
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dielectric FUNCTION(IMACINARY) 



OiELECTRiC FU«CrOR(iliWCiNARY) 



Figure 3. 7 . 9 


Imaginary part- olT the complex dielectric 
Inunction as a t^unction of* f^reguency f^or the 

model case within aCPA with c=0. 5 & 0.2 

respectively. — with correction 

+ without correction 


niElECTRii' FUNCTiQN(REAL) DIELECTRIC FUNCTIOH(REAL) 




Real part of the complex dielectric function as a 
function of frequency for the model case 

within 2CPA with c=0. 5 & 0.2 respectively. 

— with correction 
+ without correction 






REFRACTf/E INDEX REFRACTME INDEX 




Figure 3 . 7.12 Figure 3.7.13 


Refractive index function as a function of 
frequency for the model case within 2CPA with 
c=0. 5 & 0.2 respectively. 

— with correction 
+ without correction 


ABSORPTK>N COEFFICIENT 



ABSORPTION COEFFICIENT 



Figure 3. 7. 15 


Absorption coefficient func^on as ^ 

frequency for the model case within 2CPA 

with c*0.5 & 0.2 respectively. 

_ with correction 
+ wi thout cor r ec ti on 





I wo 


TRANSITION PROBABILITY 



— wuBPimnw + wiKwicmcnw 


Figure 3.7.16 


TRANSiriON PROBABIUTT 



— wti tHHErriM 't wma: Bswccnei! 


Figure 3 . 7 . 17 


Square of the optical matrix element Ctransition 
probability!!) as a function of energy for the model 
case with c= 0.5 and w = 4.62 within ICPA and 
2CPA respectively. — with correction 

+ without correction 


TRANSITION PROBABILIPf 



^ t iBiweiweiwi* 

Figure 3.7.18 


TRANSITION PROBABILITY 



-M -4 HI ^ 

Binit — > “H 

— wiBMKBtn f wummucam 

Figure 3, 7. 19 


Square of the optical matrix element Ctransition 
probabiiityD as a function of energy for the model 
case with c= 0.5 and o> = 9.24 within ICPA and 
2CPA respectively. — with correction 

+ without correction 






TRANSITION PROBABILITY 



Figure 3 . 7 . £0 


TRANSITION PROBABIUTT 



Square of the optical matrix element C transit! on 
probability!) as a function of energy for the model 
case with c= 0.2 and <0 = 4.62 within ICPA and 
2CPA respectively. — with correction 

+ without correction 


TRANSITION PROBABILITY 


iMi 



TRANSITION PROBABILITY 



Figure 3.7.23 


Square of the optical matrix 
^obabilltyj as a function of 

Lse with c= 0.2 and « = 9-24 within ICPA and 

SrPA respectively. — with correction 
2CPA respectively ^ ^^^hout correction 






OOS*TRANSITION PROftikSILITT 



Figure 3.7.24 


OOS*TRAHSiTION PROWILITY 



Figure 3. 7. 25 


Producl of the valence band DOS and square of 
the optical matrix element Ctransition 
probability^ as a function of energy for the 
model case with c= O. 5 and <»> = 4. 62 within ICPA 
and 2CPA respectively. — with correction 

+ without correction 


OOS^IRANSITION PROB^BIUTY 



DOSnRANSITION PROBABlLIPr' 



Product of the valence band DOS 

the optical tton^^^f’^ energy for the 

probability^ as a ^ = 9.24 within ICPA 

model case with c- 0.5 ano o correction 
and 2CPA respectively. -- correction 
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DOSnRANSITION PROMILITY 



Figure 3.7.28 


DOSnRAHSmON PROWILITY 



Product of the valence band EX>S and square of 
the optical matrix element C transit! on 
probability^ as a function of energy for the 
model case with c= 0.2 and <o = 4.62 within ICPA 
and 2CPA respectively. — with correction 

+ without correction 


WS*rRANSIT)ON PROBABiLiTt 



D0S*TRAN5ITI0N PROBABILITY 



Figure 3.7.31 


Product of the valence band IXS and square of 
the optical matrix element Ctransition 
probability^ as a function of energy for the 
model case with c^ 0.2 and co = 9.24 within ICPA 
and 2CPA respectively. — with correction 

+ without correction 




CHAPTER IV 


OPTICAL PROPERTIES OF III~V TERNARY ALLOYS 


4. 1 INTRODUCTION 


The application of the Augmented Space formulation of 
optical conductivity on a simple model for random-semiconducting 
binary alloys discussed in the previous chapter gave encouraging 
results and it provided us with sufficient theoretical basis for 
applying the formulation to more realistic situations of III-V 
ternary alloys. 

The III-V semiconductors GaAs, InAs, GaSb, GaP and InP, 
etc. form zinc-blende structures Csee figure 4.1. ID. A zinc-blende 
structure consists of two interpenetrating f.c.c. sublattics, the 
ariion Ce.g. As> and the cation Ce. g. GaD sublattices, which are 
displaced from each other by a vector t = Cl/'4, 1/4, 1/4D a, where a 
is the lattice constant. Consequently each anion in the crystal is 
surrounded by a tetrahedral arrangement of four cations and vice 
versa. A schematic picture of the flattened network of the 
structure is shown in figure C4.2.2D. 

Pairs of these semiconductors mix in all concentrations 
form ternary substitutional disordered alloys. The sublattice of 



I j 




Zinc blende crystol structure 

Figure 4.1.1 
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the common ion is ordered while that of the other ion is occupied 
randomly by one of the two possible constituents. 

The work reported in this chapter has two parts. In the 

first part, as we adopt basically a tight binding view point, we 

attempt with the help of the Recursion technique C§4. 3;> to develop a 

fully real space method for generating the various constituent Green 

function for GaAs, GaSb, and InSb and the alloys GaAs^_^Sb^ and 

Ga In Sb Cfor c = O.B and O. SD. We then use the CPA to generate 
*-c c 

the valenece band density of states for the above alloys. The basic 
Hamiltonian is set up in the bond orbital basis C§4.2D with 
parameters fitted to experiment. 

In the second part we apply our formulation for optical 
conductivity developed in the previous chapter to calculate various 
optical response functions for the above alloys. 

I I ^ 

The modified square of the optical matrix element, 19;^^ I 
is calculated with the use of known numerical values of wave 
functions for different atoms reported by Herman and Skillman 
C1970D. Again the wave functions are written in the bond orbital 

basi s. 


4.2 HAMILTC»(IAN IN BgKD C^BIT^ BA^S 

To illustrate what are bond basis, let us first take a 
diatomic molecule. The system is formed by 


simple example of a 
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two atoms A and B Csee the figure 4.2.15. 


The LCAO associated with A atom is - r^5 and that of 

with B atom is — r g5 . Then the secular determinant in the 

LCAO basis is given by 


”aA “ ^AA 


“bA “ ^ ^BA 


»AB - - ^AB 


*^BB ^ ^BB 


C4. 2. 15 


where Hamiltonian H in LCAO basis is 


H = 


”aa “ab 


”ba ^BB 


C4. 2. 25 


”aA ' " ’'a^ H .^r - r^5 dr = 


Hgg » S 4,*Cr - rg5 H ^ Cr - rg5 dr = 
«BA “ «AB * ^ 

^AA “ ” ^A^ <f> Cr - r^5 dr = 


C4.2. 35 


1 



BB 


x <p cF - 


^ Cr - TgZ) dF = 1 
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^AB “ ^BA 


X <P Cr ~ r^:) <p Cr - F^D dF 0 


energies, V is interaction potential 
between A and B atoms and ^ is the eigenvalues of H. Therefore a 
representation of H in the LCAO basis is now given by 


H 


V 

£ 


B- 


C4. £. 4) 


As an example, the solution of the secular equation for a simple 
case where “ '^B “ are and where 


s » 


£ . * 


£ ~ V 

1 - s 

c + V 

1 + s 


C 4 . £. 5:5 


The* correspondi ng eigenfunctions are given by 




^ Cr5 » 


y 2C1-SD 




*s 


C4. 2. 62) 


-/eTcT^ 
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4>_Cr'> is called anti bonding state and 4>^Cr^ the bonding state. 

The representation of the Hamiltonian in the new basis C4>_, is 

diagonal, with s_ and in the diagonal positions 

The work by Harrison and his coworkers Cl 973-74:) has 
demonstrated that the systematics of bonding of various physical 
properties can be described by the bond orbital model CBOMZ) in 
covalently bonded systems. In semiconducting materials the set of 
bonding states spans the valence band and the set of anti bonding 
states spans the conduction band, if the hybridising term between 
them is negligible. Thus the advantage of working in this basis is 
that if we want only the valence band structure, we can only work 
with the half of the space spanned by the bonding states. These 

3 

basis functions are constructed by linear combinations of sp 
hybridised orbitals on adjacent atoms. 

To get the Green operator corresponding to bonding states 
C valence band!!), we may partition the Hamiltonian as 




a ^ ii - - H 


C4.a. 73 


wh^r^ And proj€fC't on'to ■the vAlence Cor bondinQD and conduction 
Cor anti bond! ng> subspaces respectively and is the hybridising 

term. 


gyCzD « t*L “ liv " - cv ^ 


C4. 2. 8D 
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G CzD = [ zl. - H ] ~^c 


C4. 8. 91) 


The valence band density of states is given by 


- Ciyn^ Ta, 


C4. a. io:3 


Within the bond orbital model, it is assumed that H is 

— cv 

small , and so we may completely neglect the antibonding part and 
work exclusively in the bonding subspace to get valence band density 
of states. In case there is significant coupling between bonding 
and anti bonding states, the bond orbital approximation is 
inaccurate, and the correction terms have to be included. 


Using the notations of Chen and Sher C1978D we fix the 
anion f.c.c. sublattice as the reference lattice on which each 
lattice point is represented by a lattice vector J. The four 
bonding orbitals surrounding it and pointing in the tetrahedral 
directions, may be represented by |ja>, a = 1,£,3,4 centred at 

points specified by J + where = Cl, 1,15 a/8, Tg = C-1,-1,15 
a/0, a C-1, 1, -15 a/8 and r = Cl, -1,-15 a/8, where a is the 

lattice constant. The bond orbitals themselves are composed of 
symmetric linear combinations of the hybridised atomic orbitals |a> 
and jc> centred at the anion and cation atomic sites respectively, 
i * 


> » Cl/N^> t |a> + |c>] = |ac> 


C4. 2 . 115 
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WllSF €? N is "Lfl^ nOFlTLal i 4 ^ ± . j — 

« normalization constant y^2Cl+SD CS = <a|c>, the 
orbital overlap of anion-cation pairD, and 

|a> = ci/a> [s^> + cy5/2:> |p^> 

C4. 2. 12:) 

[c> = Cl/25 (s^> + c-/3yz^ |p^> 

where {s^> and |p^> are the s and p orbitals on the anions while 
[^c** and are those on the cations. The p-orbitals are taken 

to be directed along the bond in the above expressions. When 
expanded in terms of a bond basis set <|J,o(>> the BOM Hamiltonian 
1 s gi ven by 

= E E E E |J»«> <J.C([H^ jJ'ot'Xj'a' { C4.2.13:) 


If we consider bond to bond overlaps upto second nearest neighbours 
bonds, there are five distinct overlaps <J,o{ which are 

defined below Csee figure 4.1.23. 

D : The bond energy <j«|H 

: The matrix elements between adjacent bonds with an anion 

in common. 

y® : The matrix elements between adjacent bonds with a cation 
1 

in common. 

D : The matrix elements between parallel second nearest 

z 

nei ghbour bonds . 
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: The matrix elements between non-parallel second nearest 

nei ghbour bonds . 

Chen and Sher C1978D have shown that the inclusion of 
thi rd-nearest neighbour interactions, does not change the band 
structure significantly. So, retaining the interactions upto second 

nearest neighbour parallel and non— parallel bonds only, the distinct 

matric elements are D, , D and y 

1 1 211 2 

4-3 THE RECURSI METHOD 


Band theory for systems, like simple metals where the 
interaction between electrons and atoms is weak, exploits the 
physics of long range periodicity of the atomic potentials to 
express electronic properties of the solid as a coherent 
superposition of the electronic properties of all the atoms. When a 
strong interaction is present between the electrons and the atoms, 
the electronic properties no longer depend on long range periodicity 
and thus this picture of band theory breaks down. In such cases the 
properties depend only on the first few shells of neighbours of each 
atom. The d electrons in transition metals are good examples. While 
the band theory is still a valid formal solution to the Schrddinger 
equation, the physics is better understood by means of a solution 
that explicitly accounts for the role of the local environment. 


To do this, we have to calculate the local density of 
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states. The local density of states describes the effect of the 
rest of a solid on one region. Physically, the local density of 
states is the intensity of each eigenstate on a particular atom or 
bond. Mathematically, it is the magnitude squared of the projection 
of each eigenstate on a local orbital. Most properties of solids 
including binding energy, photoemission spectra, and magnetism can 
be related to the local density of states. 

The basic problem is to find the local density of states 
from a ti ght -binding or localized orbital model. The physics 
demands that it be done in such a way that the local orbital itself 
has the greatest effect, and that successively more distant orbitals 
have lesser effect. The solution must define a hierarchy of 
environments so that their relative influences areexplicitely 
displayed in the local density of states. 

The Recursion method, introduced first by Haydock 
ei'.ai.Cl973> and discussed in detail by Haydock C1980D, provides one 
of the best real space methods to tackle the above problem. In this 
method, linear combinations of the local orbitals are so determined 
that the electron passes through each state as it diffuses away from 
the local orbital from which it starts. Thus it serves the demand 
of the physics Just described above. 

The sequence of environments are described by orthogonal 
slates IV orbitals, 

lb practice tha states |u„> are localised on the shell of atoms n 



laa 


hops from the atom accommodating [<^^>. if the local density of 
states is to be obtained around a particular atom or bond described 
by a basis orbital > |u^> is taken equal to |^^>. the rest of 

the orthogonal states are generated from the condition 




+ b 


n+1 




u . > 
n -^1 


C4. 3. ID 


where H is the Hamiltonian of the model. a^ and b^ are the coupling 


constants of each environment to Itself and its neighbours. 


The coefficients and b^ are generated recursively from 


u> = {d'>i. b =1 

I o I o 


» 

n 


<u Hfu > 
n — n 


<. u u > 
n ' n 


C4. 3. ED 


n+l 


b <u u > 
n n * n 

<u fii' > 

n+l n+1 


We have used the orthogonality of the new basis <|Uj^>>. 
The above equations are a convenient algorithm for computational 
generation of the coefficients a^ and b^. It is clear, from the 
equation C4.3.1D. that H has a tridiagonal representation in < lu^» 


basis : 
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Equivalently, given the Hamiltonian matrix H in orbital basis 

<|^^>>, the Recursion method generates a unitary transformation U 
such that 


U H y “ C4. 3. 42 ) 

where s tri diagonal. 

In a semi -infinite chain model described by the orbitals 

<(u >> and the parameters a and b .the Hamiltonian has a 
' n n n 

tridiagonal representation and is symmetric in < [u^>> basis CHaydock 

1 980 > . 


Therefore, in other words the Recursion method transforms 
a complicated three dimensional model into semiinfinite chain model 
described by the parameters a^ and b^, thus making the problem much 
more simpler. A graphical representation of the chain model is 
shown in the figure C4.3. 13. 

Let us now get the expression for the local density of 
states n CE>. The corresponding Green function is given by 



hond 


origin 

o , 

A c 



Fi gur 4. a. 1 


Co ordinates of molecules in a diatomic 
molecule with respect to the origin atT^S 
centre of the bond 





t « # 


l«3> 


A graphic representation of the chain model 
including the states 

recurrence parameters <a > and <b > 

n n 


Fi gure 4. 3. 1 
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G^Cz^ = <<PjCzl - K) ^ \<P > 


<u (czL - ro ^ j u > 
o' — — ' o 


<u 


Cz-a!) 

o 

-b 

1 

O 


-b 

1 

Cz-a D 

1 


-b Cz-a D 
2 2 


o o 

b_ O 

b 


u > 

o 


C4. 3. 5D 


Let us define the determinant of the matrix with the first 

n rows and n columns as D CzD, then 

n 


G Cz:) = [D CzD/ D CzD] 
o 1 o 

= CD CzD/<Cz-a 3D Cz:) -b^DCz3>] 

i O 1 12 


= l/CCE-a 3-b ^ D Cz3/D Cz3 ] 
0 12 1 

= l/CCE-a 3 -b ^ G Cz3 3 
o 11 


If we continue this process upto n-steps^ we have a 


continued fraction expansion for the Green function 
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z - a 

o 


1 



b* G CzD 

n n 


where 

e Cz5 = l/[z-a -b^ G CzD3 C4.3.7D 

n n n+1 n+ 1 

The continued fraction representation, in general, may not 
terminate after a finite number of steps. It has been verified that 
for many systems the continued fraction converges very fast. When 
the convergence is slow or oscillatory it can be terminated by 
various termination schemes. Several termination schemes have been 
developed by Nex Cl 978, 1984D and Haydock and Nex C1984Z). 

The density of states n^CED is now given by 


n CE3 = - Cl/nD Um, ,S‘’iTv <3 CzD 

o 71-»g410 o 


= S'4n. 


-Cl/nD 



E - a - b^ G CE+iOD 
n n n 
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4. 4 THEORY OF C(>fTlNUED FRACTION TERMINATORS 


To get a proper terminator scheme to truncate the infinite 

continued fraction expansion, one needs first to examine the 

asymptotic behaviour of the continued fraction expansion. In the 

continued fraction expansionC4. 3. 73 , one examines relations between 

features of the local density of states n CE3 and the coeffi cents 

o 

<a .b >. 
n n 

The asymptotic description of a^ and b^ is based on 
two elements : 

Ci3 the main asymptotic behaviour depends on the band structure 
, i. e. ^^n^ converges towards limits in the single 

band case or oscillate in a predictable way in the multi band 
case. 

Cil3 Damped oscillations are created by isolated si ngul atari ties 
of n^. The period of oscillation is related to the position of 
the singularity, and the rate of damping is related to the 
nature of the singularity. 

The relation C4. 3. 13 defines |Uj^> as a polynomial in H 
acting on |u^> : 

IV “ Pn IV C4.4.13 



I 


where, 

p CH5 = CH - a 5/b etc. C4. 4. a:5 

1 ~ o i 

These polynomials depend only on the sequences ^^n^’ they 

are generated by the recurrence relation 


b p CzD = z-a 

1*^1 o 

C4. 4. 3:> 

b p CzD = Cz-a D p CzD -bp CzD, n > 1 
n-M n+1 n n n n-i 

We shall not go into the detail of the behaviour of the 
polynomials p^CzD. The details have been given elsewhere C Magnus, 
1984D. Let us only concentrate on the terminator schemes. We 
denote the solution of the recurrence relation C4.4.3D by q^CzZ) with 

Cz-a :> q CzD - b q CzD = 1 C4. 4. 4D 

O 11 

Then one has 

q CzD = Cz-a - b q CzD/q CzD) ^ 
o olio 

C4. 4. 53 




-1 


. S^CzJ = ■ '’n+1 


-1 


or 


C4. 4. 6D 
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G = [p^^Vp ] + 1/Cb p Cp /G -p D] C4.4.103 

O n-1 n n n n n n-i 

Thus showing that the usual continued fraction approximant p^^^/p 

n—l r» 

corresponds to the crudest estimate G =0 for the terminator. Much 

n 

better results will be obtained with terminators coming from model 

continued fractions with a non-vanishing imaginary part. The models 

are chosen according to the behaviour of the local density of 

states, n CEZ>. 
o 

THE SQUARE ROOT TERMI NATOR 


In the case where the system has multiple bands Cin 
particular single bandZ) and n^CED has only square-root singularities 
at the band edges, the widely used terminator is The Square Root 
Terminator . Let us consider the continued fraction expansion of a 
root of the quadratic equation 

RG^-SG + V = O C4. 4.115 

o o 

where R, S aqnd V are real pol ynoml nal s . 

If the root is written as : 


G Cz5 
o 


[SCz5 - X*^*Cz5] / [aRCz5] 


C4. 4. 135 
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XCz5 = S*Cz5 - 4 RCz3 VCz5 C4. 4. 135 

= Cz-E 5CZ-E 5. . . Cz-E 5 

1 2 2m 

E < E < E also assumed to be real, G will have a 

1 2 2m o 

non-vanishing imaginary part on the region where XCz5 < O, which 
is the set made of m intervals. E and E are the band edges 

2m~l 2m 

of the m^^ band. 

The interesting feature of the form C4. 4.135 is its 
reproducibility under the transformation C4.4.65 ! 


if - 4R V = X 

n n n 


C4. 4. 145 


G = CS - X^'^^]/3R = 3 V /[S + C4. 4. 155 

n n n n n 


^ E2Cz~a DV /b 3 -S 

Z “a -b n n n n 

n n~l— 


X 


1/2 


-1 


3b V /b 

n+ 1 n n 


SO that 


R Cz5 / b 

n+i 


n-M 


V Cz5 / b 

n r> 


S Cz5 = 3Cz-a 5 V Cz5/b -S Cz5 

r» n n n 


b V Cz5 
n+i n-ti 


b 

n-»*i 


S 


Cz5/4R 


ri+ 1 


Cz5] 


C4. 4. 165 
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= Cz-a 3CS CzD - S CzD5 / 2+b R CzD 

n n-M n n n 


Moreover, the degrees of R , S and V will decrease and stabilise 

n n n 

ultimately on the values m+1 , m and m-1 . Indeed by putting the 
equation C4.4.155 in the equation C4. 4. 7:>, one obtains R and S in 

n n 

terms ofP^, after transformation using the equation 

C4. 4. q :> 


R -<b 

n n 


, 1/2 




Rq 


n-l 


C4. 4. 17:5 


S 

n 


X^'^^b Cp 

n n-± 


q D 


+ SR b q q 


As q CzD behaves likeb.,..b when z oo> the degrees of 

n n n 

R and S will decrease by steps of 2 units^ down m-1 and m. 
n n ^ 

The easiest case occurs when m = 1 Csingle band3 , where 


one has 
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R Cz3 b = 1 

n n 

a =CE+ED/a C44185 

rt-i 12 N. ». **. j. 

S Cz5 = z-CE + E5/a 
n 12 

= CE^ + 4 if n is large enough 

The corresponding orthogonal polynomials are known as the 

Berstein-Szego polynomials CChihara, 1978D. Practically CE +E D/2 

1 2 

and CE j-E D/4 must be the limits of a and b when n -» co, a very 
2 1 n n 

well known situation. 

4. 5 DENSITY OF STATES OF III-V TERNARY ALLOYS 

In this section the ICPA together with the Recursion 
method will be applied to calculate the valence band density of 
states of III-V random ternary alloys . As before the conduction 
band DOS will be taken to be like that for free electrons Csee 
figure C4. 5. ID D. 

Before calculating EXDS for alloys » the valence band DOS 
of pure semiconductors GaAs , GaSb and InSb will be calculated using 
the Recursion method . The Hamiltonians for such systems in the 
moleculer bond basis are given by the equation C4. 2.13D . The 
parameters like D . y® , D^and y^ Cdefined in § 4. 2D for such 
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systems are taken from the work by Chen and Sher C1978D . These are 
shown in the table C4. 5. ID 

Table 4. 5. 1 




Interaction 

Parameters 

CeVD 



D 

Y 

* IS 

Y 

* ±A 

D 

2 

Y 
^ 2 

GaAs 

-11.1687 

-1 . 51 87 

-0. 4350 

0. 3406 

-0. 1094 

GaSb 

-9. 6437 

-1 . 3687 

-0. 3750 

0. 3906 

-O. 0594 

GaP 

-10. 8135 

-1 . 3635 

-0. 3350 

0. 3813 

-O. 1188 

InAs 

-1 0. 0563 

-1 . 3937 

-0. 5135 

0. 3844 

-0. 0656 

InSb 

-9. 5635 

-1 . 3500 

-0. 3875 

0. 3937 

-0. 0563 

InP 

-10. 3750 

-1 . 3135 

-0. 4135 

0. 3635 

-O. 0375 


where 



and 


1 A 



Y 


c 

1 




The band gaps and the effective electron masses for the 
semiconductors GaAs » Gasb » and InSb are given in the table C4,S.2]) 
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TABLE 4. 5. a 


Material 

Band gap 

M = Cm /m 

< 


E CeVD 



9 


GaAs 

1 . 53 

0. 068 

GaSb 

0.78 

O. 050 

InSb 

O. 23 

0. 014 


Recursion parameters Ca's and b’s , defined in § 4. 3 3) for 
these semiconductors are calculated upto 14 steps C shown in the 
teble C4.5.3> D and the continued fraction expansion is terminated 
using the square root terminator described in § 4. 4 . 



136 


TABLE 4.5.3 
Recursion Parameters 


GaAs 

GaSb 

InSb 

CeV3 a 

b 

a 

b 

a 

b 

-11.16870 

15. 76212 

-9. 64370 

11. 45910 

-9. 56250 

12. 38153 

-13. 95030 

8. 29524 

-1 1 . 54485 

5. 96766 

-1 1 . 97702 

5. 82638 

-1 1 . 81 261 

11.17063 

-9. 49665 

6. 23608 

-10. 09841 

7. 00559 

-13. 10205 

17. 32363 

-9. 37953 

8. 93029 

-10. 82950 

13. 72379 

-12. 59923 

6. 40295 

-11.73653 

8. 57858 

-11. 23511 

4. 23485 

-1 1 . 63092 

18. 17576 

-8. 64858 

3. 43360 

-9. 30902 

10. 67344 

-14. 01978 

7. 27838 

-10. 12119 

9. 66969 

-12. 33525 

7. 46476 

-11 . 62677 

1 1 . 47339 

-11.13626 

8. 54560 

-9. 93838 

5. 47466 

-13. 65308 

15. 04109 

-9. 49544 

4. 30477 

-10. 85169 

12. 51796 

-12. 18002 

7. 72794 

-9. 681 94 

8. 71177 

-11. 28667 

6. 06868 

-12. 781 03 

14. 64544 

-11. 26134 

5. 63407 

-10. 19797 

8. 54836 

-12. 93502 

7. 77484 

-10. 33605 

6. 66237 

-11.39538 

7. 75767 

-12. 67082 

7. 90878 

-lO. 05651 

5. 59190 

-11. 21823 

6. 87441 

-11.85587 

10. 77620 

-11.02555 

7. 40693 

-10. 03531 

4. 78794 

The 

val ence 

band DOS per 

bond f or 

pure GaAs 

, GaSb and 


InSb are shown in figures 04.5.33 , C4.5.33 and C4. 5. 43 respectively 
.It is clear from these figures that all essential features with a 
S”*like band separated from a p—like band , which itself is split by 
the cubic crystal field into a e^ and a t^ like points . The heights 




137 


are ir> the ratio 1:1:2 as expected . There are still many 
oscillations which may be artifacts of the terminator , particularly 
in the s-like part and the upper edge of the p-like part of the band 
C V. K Srivastava , 1982 2 > . 

Now the valence band IX)S per bond will be calculated for 

the random semiconducting alloys . Two alloys will be chosen : 

GaAs Sb and Ga In Sb for c=0. 5 and 0.2 . Table C4.5.4D shows 
i-c c i-c c 

that each of the six semiconductors have nearly the same bond 
length . Hence , it will be assumed that the interaction parameters 
corresponding to the pure materials are the same as those in the 
alloy . The vaccum state sets the same zero of the energy scale for 
all the three materials . 



In an 

alloy the question of 

interaction 

between 

two 

differ ent 

kinds 

of bonds arises . For 

exampl e in 

the case 

of 

GaAs Sb 

alloy , 

a configuration of the kind 




Ga 

1/ \ s 

/ \ 

/ \ 

As Sb 

may be a part of the alloy . It is shown below that the interaction 
bonds 1 and 2 is the average of y^'s . defined in 

-i2 * 

§C4. 25 , corresponding to pure Ga-As and Ga-Sb crystals as a 

reasonable approximation . 
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Table 4. 5. 4 


Material 

Bond length 

Bond energy 

CeV5 

Wor k f unc t i on 

CeVD 

GaAs 

2. 45 

-5. 6687 

5. 5 

InAs 

2. 61 

-4. 7562 

5. 3 

GaSb 

2. 65 

-4. 7437 

4. 9 

InSb 

2. 81 

-4. 7625 

4. 8 

GaP 

2. 36 

-5. 1125 

5. 7 

InP 

2. 54 

-4. 5750 

5. 7 


Let the bonds 1 and 2 be described by the quantum states 
> and > respectively . Then 


> = 1/N 

I Ga^> 


\<p^ > = 1/N 

iGa > 

* 2 

isb^>] 


in terms of sp® hybridised orbitals at various sites directed along 
the respective bonds and N is the normalization constant rst 1[2 

.assumed to be same for the GaAs and GaSb bonds . We are concerned 

with t.e. , 
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<<P. 


^®xaci , 
Ji i 

--12 I 


= 


Cl/N^D 


f< <Ga I + <As I > < iGa > + ISb > >1 

j_ 1 ' 1 ' —12 '2 2 J 



The first term is vanishing because of the orthogonality of two 
different hybridised orbitals at the same sites . The middle two 
terms would be identical in the unalloyed case due to the 
geometrical symmetry of the location of the orbitals involved 
Hence this part of the expression is equivalent to the average of 
these corresponding to the pure cases. The last term is small 
compared to the middle two since it involes more distant orbitals . 
The approximation therefore involves the neglect of the variation 
in this term over the constituent materials . 

Therefore for these alloys interaction parameters , r^, 

D and y are given by the average of the above parameters of 

1 f 2 2 

the respective pure semiconductors . 

The band gaps and the effective electron masses for these 
alloys are taken as the average of the respective pure 
semiconductors . These and the interaction parameters are calculated 
and given in the table C4.5. 5D 
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Table 4. 5. 5 


Mater i al 

Inter action 

Parameters CeV3 



* lA 

r D 

ic 

> y Band 

2 *2 

gap 

M=C m /m D 

e 

GaAs Sb 

O. » 0,5 

-1 . 7937 

-0. 9937 

0. 3556 ~0. 0844 

1.155 

0. 059 

GaAs Sb 

O. B 0.2 

-1 . 8837 

-0. 9937 

0. 3506 -0. 0994 

1 . 380 

O. 064 

Ga In Sb 

0.5 0.5 

-1 . 6906 

-0. 9281 

0.3712 -0.0588 

0. 505 

0. 032 

Ga I n Sb 

0.8 0.2 

-1 . 662S 

-0. 9075 

0.3712 -0.0588 

O. 670 

0. 043 


Using these average values as inputs in the Recursion 
method , we calculate the EX>S , within the virtual crystal 
appaoxi mati on CVCAD , for the above alloys . Again the recursion 
parameters are calculated upto 14 steps _of the recursion and the 
continued fraction expansion is terminated using the square root 
terminator . These OOS are then used as the input for the initial 
choice of <G> in the CPA equation C3.2.22D for our CPA calculation . 
Finally new self energy E and from this new <G> is calculated self 
consistently . 


CPA 


DOS 


Fi gur es C 4 . 5. 5) > C 4 . S . 63 » C 4. 5. 73 and C 4. 5. 83 show the 
for ®aAs^ jjSb^ g . ®^^o. 8 ^o .2 ’ 
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2 ^^ respectively .The bond energies randomness is greatest 
in the former set of alloys , but not enough to cause split or 
impurity bands . Therefore as indicated by model calculations the 
50-50 GaAsSb shows the maximum effect of randomness . 

4. 6 OPTICAL PROPERTIES OF III-V RANOOM TERNARY ALLOYS : 

RESULTS AND DISCUSSION 


Having calculated the valence band IX>S for GaAs Sb and 

1 - c c 

Ga In Sb Cfor c = O. 5 and 0. 2D ► we shall now apply the Augmented 

1 ■“ C C 

space formalism to these alloys . Let us first write the centred 
orbitals for different species in the alloy . Take the example of 
GaAs Sb . The spacies A is GaAs and B is GaSb . Therefore 

1 - c c 

\<p > = \<p > - Cl/y§D [[<;«> > + \<p >1 C4.6.1D 

i^A '^OaAs |_'^Oa I^As J 

where I*^as^ centred orbitals for the respective 

3' 

atoms and are given by the linear combinations of sp hybridised 
or bi t al s C LCAOD : 

= C1/2D + CV^/2D and 

C4. 6. 2D 

1^*.^ = ci^3 |,^•^> cys^5 

Where |4>* > and \(p^ > refer to the s and p states for Ga . Similarly 

1^* > and \4>^ > are defined . Therefore 
'^A» »^A» 



DOS“^ 


1 42 


CONDUCTION BhND DOS 



Conduction band DOS Cto 
be free electron likeD 
as a function of energy. 


mEMCE BAND DOS 



Valence band DOS for the pure 
GaSbCupto 14 steps of recurs! onD 

as a function of energy. 


WfNCE BAND DOS 



Valence band DOS for the pure 
GaAsCupto 14 steps of recursion^ 
as a function of energy. 


VALENCE BAND DOS 



-a —-a -n Hj 

-17 Ecev; -s 


Figure 4. 5. 4 


Valence band DOS for the pure 
InSbCupto 14 steps Qf recursion, 
as a function of energy. 
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Figure 4. S. 6 
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= ci-'vSj [ci/a:>[|^*^> + |^•^>] cyS'x2>[|*J^> |«p^>]] 

= ci^vs:> [[|^•^> |^■^>] ^ cyS>[|^J^> + |^'.>]] 

Si mi 1 ar 1 y 

\ 4 ,^> . ci^yg:. [[|^;„> |^•^>] + cyS3[|4iJ^> + Kt>]] 

The Fourier transforms of <p CrD and <j> dr'i are given by 

A ® 

\4> CkD> = r e’'’^'''^ CFD d®r 

* A,B ^ A,B 


= Cl /EV^D [ 

(j!>® Ck5 

+ <*® Ck5 ] + 


cy^/2yS5 [ 

4>^ CkD 

^Ga 

+ 4>'^ ^Ck5 ] 

C4. 6. 45 

where ^®Ck5 = CVn/EkD 

00 

f r sinCkr/E5 <^®Cr5 dr 



O 


C4- 6. 5:) 

00 ^ 

<p^Ck':> = Ct/k*DV§Ti J [ sinCkr/2D - Ckr/a5 cosCkr/2D ] <^^CrD dr 

o 

Cr in Bhor unitD 

k==2-|MCE^''+ o>5" 

m 

Here E and w are in Rydberg units and M = m /mg.- 

i 

Thus using the equation C3.6.323 and the calculated 

n CE 5 , n CE!> and self energy together with the equation C4.6. 35 
i i f 


DOS 
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C4.6. 43 and C4. 6. 53 and the values of <^*Cr3 and ^*'Cr3 C reported by- 
Herman and Sklllman C197033 for the various atoms in the alloys , we 
calculate the optical conductivities for the alloys . Since the 
values of ^ Cr3 and 0*^Cr3 are given in the unequal intervals of r 
Cin the logerthmic scale3 and are zero after a finite values of r , 
we used trapezial rule to calculate the integrals C4.6. 53 . Then 

using the Kramers— Krdnig relations , we calculate the other response 
functions 

Figures C4.6. 13 and C4. 6. 33 show the optical 
conductivities as a function of the incoming photon for GaAs Sb 

0.5 0.5 

and GaAs ^Sb respectively . Each of these figures show 

significant effects of disorder . It is clear that in the case of 
50?^ - 50% concentration the effect is more than the 80% - 30% 

concentration case . In both of the cases the effect of disorder has 
reduced the magnitude of the conductivity . Thus we can say that 
the vertex correction is negative . 

The calculations of the other response functions like the 
imaginary part of the complex dielectric function Cfigures C4. 6. 33 
and C4. 6, 433 .the real part of the complex dielectric function 
Cfigures C4.6. 53 and C4.6. 633 .the refractive index function Cfigure 
C4. 6. 733 and the absorption coefficient Cfigure C4.6.833 also show 
the significant effects of the disorder .These quantities are 
related to one another as discussed . In these figures we see that 
the sharp structure of the DOS gets smoothened out because of the 
fact that in these properties DOS is con-snjlated with featureless EX>S 
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of conduction band and relatively featureless of the square of the 
optical matrix element .These figures also show that peak positions 
are shifted and widths are increased due to the disorder induced 
corrections . 

Figures C4. 6. 9D . C4.6. lO!?, C4. 6. 115 and C4.6.ia5 show the 
square of the optical matrix element Ctransition probabilityD as a 
function of energy for GaAs Sb Cfor = 0. 3 and 0. 9Ry5 

O • 9 O . is 

and GaAs _Sb Cfor w = O. 3 and O. 9Ry5 respectively .It is clear 
from all these figures that optical matrix element indeed depends on 
both the energy and disorder . Though in the case of C80-20%5 
composition, the magnetude due to disorder effect does not differ 
much CThe effect is least in the weak disorder case of the dilute 
alloyD but there is some sharp structures seen in this case. These 
structures are more in the case when a> = 0. 9 Ry. This is because 
the impurity CGaSbD has a site energy 2 = 0. 7 Ry which lies between O 
to 0.9 Ry and the scattering is more at impurity sites. Thus this 
clearly shows that in order to interpret energy resolved 
photo-emission experiment one must take into account of the effect 
of disorder, because the energy resolved experiment gives the EXDS 
multiplied by square of the optical matrix element. Therefore the 
observation of photo— emi ssi on experiment will differ from the 
correct band structure by both the magnetude and structure. 

The figures C4. 6. 135 to C4.6.165 for G&As^ ^ Co = O. 3 

and 0.9 Ry5 and GaAs Sb Co = 0.3 and 0.9 Ry5 respectively show 

O « tEI 0*2 

how the effect of disorder make the difference of band structure 
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Figure 4.6.15 
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Figure 4 . 6 . 16 
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calculation from that of obtained from experiments. Thus to get 
accurate IX>S one has to divide the results obtained from the 
experiments by the correct value of square of the optical matrix 
el ement . 


Figures C4.6. 17D and C4.6. 18D show the optical 
conductivity and the dielectric function CrealD respectively for 
Ga In Sb. In this case the vertex correction Is almost 

0.5 0.5 

negligible. This is because of the negligible disorder effects, as 
can be seen from the table C4.5. 15 that various parameters of GaSb 
andlnSb are almost the same, means the alloys of these at various 
concentration make a weak scattering case, i.e. effects of disorder 
over the VGA is very small. The figures C4.6. 215 to C4. 6. 245 for 
self energy as a function of energy for GaAs Sb and 

Ga^ 5 ^ further make the above reasoning clear. The self energy 

for Ga In Sb is very-very smaller than the case of GaAs Sb 
Thus the effects of disorder in the case of alloy GalnSb is 
negligible. 

Since the correction term to the optical matrix element is 

CZ /<5£5<5<* and if Z , de and 64 > are very-very small, this term has to 
o o 

be neglected which is true in the case of the alloy GalnSb. 
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Figure 4 . 6 . 20 
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Self energy CReal and Imaginary parts 
respecti velyD for GaAs Sb as a function of 

0 . 50-5 


energy. 


SELF ENERGY (REAL) 5ELP (|mac.INARY) 




Self energy CReal and Imaginary parts 
respect! velyD for 5^^o 5^ ^ function of 


energy. 





CHAPTER V 


CONCLUDING REMARKS 

In t^his work we have carried out. the application of the 

Augmented Space Formalism CASFD to a formulation of the optical 

response of random semiconducting alloys. We have shown, through 

the application of the formulation to a model system and III-V 

ternary alloys GaAs Sb and Ga In Sb Cfor c = 0. 5 and 0.23, that 

i—c c 1 — c c 

the modifications due to the effect of disorder may be significant. 
This effect comes from the inclusion of the vertex correction 
related to the averaging of the square of the optical matrix element 
together with the valence band DOS, £.e. < |g. , 1^ n.CE3> = 

<n^CE^3>. Therefore the assumption which was usually made in the 
earlier works that the optical matrix element is independent on 
disorder or very weakly dependent may not be correct in several 
situations. Thus we have been able to extend the ASF in 

generalizaing the study of optical properties of random 

semoconduct i ng al 1 oys . 

In this work we have also included cluster coherent 
potential approximation Conly in the model calculation^ and shown 
that cluster effect is indeed seen in the optical matrix element. 

In this work we assumed the conduction band density of 
states to be like nearly free electron. This needs modification. 
The conduction band should be treated more accurately than the 
nearly free electron model. Within the tight -binding model for III-V 
semiconductors this has been attempted by Chen and Sher C19803. 
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We have also Ignored the vertex correction related to the 
averaging of the conduction band DOS together with the valence band 
DOS and the square of the optical matrix element, i.&. 


<|g'”|*s 6 

' I f y ' V c 


> = 


G > |g*”|*< G > < G >, 

' ify' V -c ' vfy' V c 


Thus the vertex correction term 


<GG>=<G><G> +<G>A<G> 

VC V C V — c 

has to be included maintaining the word identity A = <5Z/<5<G>. 
This can be implemented in the Augmented Space Formalism and has 
been discussed in some detail by Thakur et al . C1987D and Mookerjee 

C1986D. 

In our calculations of the EXS we went upto the 14*’^ steps 
of recursion and used the square root terminator to terminate the 
continued fraction expansion. But this is insufficient for more 
accurate calculation. Therefore the Recursion method has to be done 
more accurately, specially correct terminator scheme has to be used. 
This has been done lately by Haydock and Nex Cl 984, 1985D. 

In our calculation we used the band gap and the effective 
electron mass for the ternary alloys to be equal to the average of 
these of the respective pure semiconductors in the alloy. But this 
approximation is not good. One has to take the correct band gap and 
effective electron mass . It has been found empirically by A. G. 
Thompson and J.C. Woolley Cl 9675 that the band gap E^Cx5 of a 
ternary alloy varies with concentration x as follows : 
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E C x2> = E + bx + cx* 

9 gi 

where E^^ is -the band gap of the lower band gap binary and b and c 

are constants withE =E +b + c; E is the band gap of the 

gz gl ’ g2 ^ 

higher band gap binary. For the alloy, GaAs Sb ,the constants b 

X X“1 

and c found, at T = 300 °K, to be equal to -0.32 eV and 1.005 eV 
respectively by G. A. Antypas and L. W. James C1970J. 

We used the values of the interaction parameters for the 
pure semiconductors obtained by parametrization procedure C which is 
very ambiguous!) of Chen and Sher. For more accurate calculation the 
value of these parameters obtained by the Chemi cal -pseudopotenti al 
approach of Anderson C1969D and Bui let t C1975J may be used. 

In our calculations we adopted tight binding model. We 
should, however, not stress too much on this, as the tight binding 
model calculation should indicate trends rather than really 
quantitatively reliable estimate. Therefore the need was felt, for 
future work, to go beyond the tight binding approximation. For 
metallic alloys, the self-consistent KKR-CPA techniques have been 
developed. The KKR equations have close resemblance to the tight 
binding method. the diagonal term in the tight binding 

Hasmiltonian is replaced by the inverse t-matrix, the off-diagonal 
hopping term V_ is replaced by the structure function , and the 

Green function by the path operator. 
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One has to go beyond the single site approximation and 
consider the cluster approximation for the correct study of band 
structure. The generalization of the Augmented Space CCPA to KKR 
method has recently been carried out by Hooker Jee C1987D. 

In cases, where the basis is not countable, an embedding 
technique was developed by Hooker jee and Bardhan 01989!) generalizing 
the earlier work by Inglesfield C1981D. It is essentially a 
generalization of the partition theorem. The Augmented Space method 
then allows us to develop a CCPA procedure in such cases. This 
work is in progress and will enable us to study systems with random 
distribution of extended impurities or defects. 

CXir calculations are at T = O °K. However, most 
experimental observations are not at T = O So it could be 

interesting to study the optical response with varying temperatures. 

Since the lattice parameters at T =300*^K for GaAs and GaSb 
are very different from one another C5. 65 A*^ and 6.10 , there is 

a lattice mismatch in the alloy of these. In other words the alloy 
of these semiconductors will have high lattice deformation at the 
boundary between these. Because of the high lattice deformation, 
these alloys are not preferred for making devices and semiconductor 
lasers. Thus this alloy was not popular, and therefore, probably, 
very little experimental studies have been carried out. Hence at 
present, our calculation can not be compared with experimental 
results. But these alloys were grown up and phase diagrams were 


163 


studied extensively CKressel and Butler Cl 9775:). Now there is 
growing interest in this alloy because it can be used to produce 
lattice matched hetero junctions in conjunction with AlGaAsSb and 
also this can be prepared most easily by Liquid Phase Epitaxy CLPE5 
method CKressel and Butler Cl 9775 5. Because of the large lattice 
parameter change with composition C7. 5% across the full composition 
range5 , a miscibility gap may exist for intermediate composition 
values. However, excellent devices have been prepared in the 1.1 eV 
bandgap energy range in the As-rich composition large CKressel and 
Butler Cl 9775 5. 

The GaSb and InSb are continuously miscible , but there 
are technological difficulties in producing homogeneous Ga^ ^In^Sb 
alloy with a given mixing ratio x by procedure which rely on 
diffusion process. The alloy has been prepared and the optical 
properties were investigated in number of studies C Madel ung , 1 9645 . 

Since in this work we are studying the disorder effect of optical 
matrix element, we shall not be much interested in this alloy , 
because the vertex correction arising due to the disorder dependence 
of the optical matrix element is negligible in this case. 
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